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Abstract 

This paper concerns the validity of the Prandtl boundary layer theory for steady, 
incompressible Navier-Stokes flows over a rotating disk. We prove that the Navier 
Stokes flows can be decomposed into Euler and Prandtl flows in the inviscid limit. In so 
doing, we develop a new set of function spaces and prove several embedding theorems 
which capture the interaction between the Prandtl scaling and the geometry of our 
domain. 
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1 Introduction 


We consider the steady incompressible Navier-Stokes equations on the domain 0 = (0, 0o) 
{Ro,oo) in polar coordinates. The boundary dQ then consists of three components: {w = 
6*0}, {td = 0}, {r = Rq}. In cartesian coordinates, the equations read: 

UU^ + VUy + Pcc = eAi/l 

UV^ + VVy + Py = eAI/ I in 0 (1.1) 

Ux -\- Vy = Q. j 

In polar coordinates, the equations read [KCOdt Page 739]: 


uu^ , , uv ^ P^ , Ur , 

-h VUr H-^- = e[Urr^ -^-^ 

ry > \ ry >^ 

C/K 


[/2 /II/ 

-^+VVr - +Pr = € Vrr + -Vr+ '' 
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fy>Z 


9 9 

j»Z 


Uu, + dr{rV) = 0. 


in Q 


( 1 . 2 ) 


Here, U and V represent the horizontal and vertical velocities of the flow, and U, V 
represent the angular and radial velocities of the flow. The Navier-Stokes equations are 
taken together with the no-slip boundary conditions on the boundary {r = Rq}. We suppose 
that the disk of radius Rq is rotating counter-clockwise with a constant angular velocity of 
Ub > 0. The no slip boundary condition in our case is then 

U\r=Ro = Ub and V\r=Ra = 0. (1.3) 

The boundary conditions at {oj = 0} and {u = 0o} will be prescribed in the text. We study 
the limit as e —)> 0. Formally, one expects solutions to the above Navier-Stokes equations to 
converge to solutions of the Euler equations with e = 0, but this does not happen due the 
mismatch at the boundary between the no slip condition enforced for solutions to Navier 
Stokes equations and the no normal flow condition enforced for solutions to Euler equations. 

To account for this mismatch, in 1904 Ludwig Prandtl proposed the formation of a 
boundary layer of size ^/e near the boundary, such that the Navier-Stokes flow can be 
decomposed into the sum of the Euler flow and the boundary layer flow. This is regarded 
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as one of the most important ideas in fluid mechanics in the last century, and the theory 
has led to astounding developments in the applied sciences. Indeed, many phenomena in 
fluids such as wake flows and plane jet flows are described by the Prandtl theory |Schflfl] . 
Despite this, a rigorous mathematical justification of the boundary layer theory remains 
open in general. 

For unsteady flows, there are several interesting results, see for instance |SC98) . |SC98) . 
|Mael4j . |As91j . |MT08j . The work of Guo and Nguyen, |GN14j . is the first result establish¬ 
ing validity of the boundary layer expansion for steady state flows in a rectangular domain 
over a moving plate. They do so using a combination of energy estimates, elliptic estimates, 
and a new positivity estimate obtained via the vorticity multiplier dy ~ (i^)' 

The main goal of this paper is to generalize Guo and Nguyen’s method in the presence of 
geometric curvature effects in order to establish the validity of the boundary layer theory 
for steady flows over a rotating disk. 


1.1 Boundary Layer Expansion 

We denote by to be an outer Euler shear flow which is radial: 

(1-4) 

Such a shear flow describes an Euler fluid which rotates counterclockwise. On the boundary 
{r = Rq}, we denote by = u^{Ro) and assume that Ue > 0. We also suppose that the 
disk of radius Rq is rotating at an angular velocity Ub > 0. We now scale to boundary layer 
variables in the following way: 

Boundary Layer Scaling: Euler Scaling: 

R = R(^r) = Rq + ^, r = r{R) = Ro + ^/i{R-Ro). 


Note that r, R > Rq > 0 and that djir{R) = yT, drR{r) = Ijy/e- We scale to boundary 
layer velocities and pressure in the following way: 

U\oj,R) = U{uj,r), V^{cj,R) = ^V{uj,r), R) = P{u,r). (1.5) 

The boundary layer velocities and pressure satisfy the following scaled Navier-Stokes 
equations: 


r 




r-U% 

+ Ve— + ' 

r 






( 1 . 6 ) 
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, ...... 1 , 1.. ... , rVR , n rU^ 

+ dR{rV^) = 0. 

We start with the following formal expansion: 

C/^w, R) = u^{io, r) + v^{(jJ, R) + y/€ul{uj, r) + \/eUp(u;, R) + e'>'’''2ri7rf, i?), 

V^{ijj,R) = Vp{uj,R) +vl{w,r) + ^/eVp{uJ,R) + /?), 

P7a;,i?) = P0(r) + P°(a;,i?) + x7Pi(a;,r) + VeP^{oo,R) + eP^ioj, R) + P^iuj, R). 

(1.9) 


(1.7) 

( 1 . 8 ) 


According to the expansions (1.7) - (1.8), the Prandtl decomposition up to leading order 
is then: 


U{uj, r) = ^{lo, R) ^ Ug(a;, r) + Up(a;, R), 

V{oj,r) = y/eV’^{oj,R) « y/(iVp{oj,R) + y/ev].{uj,r). 


( 1 . 10 ) 

( 1 . 11 ) 


[Up,ul,Up] and [Vp,v\,Vp] are approximate boundary layers to be constructed, after which 
the remainders must be constructed and controlled. We insert the expansions 

into the scaled Navier-Stokes equations, and obtain the different orders of the errors 
and K", which are detailed in equations (3.1) - (3.19). The scaled divergence free 


condition is enforced at each stage of the expansion. So, for example, for the Prandtl-0 
layer, we enforce Up^ = —dR{rVp) = —^/eVp — rUp^, and for the Euler-1 layer, we enforce 
= —dr{rv\) = — Ug — rUg^. We define the following notation which will be in use 
throughout the paper: 

Definition 1.1. Us = u^e + '^p + \7^^7 Vs = Vp + ^^Ug, and Uapp = Ug + Up, Vapp = Vs + v^. 

Once the velocities of each layer has been constructed, the pressures are defined using the 
radial error contributions up to and including the e® contributions. The e~^ order equation. 


(3.7), for instance, dictates that the initial Prandtl pressure is constant in R. The e 


order error, (3.8), is the Euler-0 pressure as given by the Euler equation for the shear radial 
flow Ug(r). We then estimate the error caused by this definition in the angular equations. 


1.2 Boundary Data 

The no-slip boundary conditions at {r = Rq} must be enforced for each order of the 


expansion in (1.7 1.8). Since the outer Euler flow is given, we have: 
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Boundary Conditions on {r = i?o}: 

uI{Rq)+ Up{u:,Rq) = Uh, Ug(a;,i?o)+ np(w,^o) = 0, i2o) = 0, (1.12) 

vl{u,Ro)+ vI{uj,Rq) = Up(u;,i?o) = 0, u^(a;,/?o) = 0. (1.13) 

Boundary Conditions on {w = 0}: 

Up(S),R) = uo{R), ul{Q,R) = ui{R), ttg(0,r) = U 6 (r), (1.14) 

vli0,r) = Vboir), u''{0,R) = v’'{0,R) = 0. (1.15) 

Boundary Conditions on {w = 0o}: 

'uli^o,r) = Vbi{r), (1.16) 

= 0 and P^r = 2 eri^. (1-17) 

Boundary Conditions as r —)■ oo: 

Up{uj, R),Up{LJ, R), [ui, vi](u}, r) —> 0 as r, i? —^ oo. (1-18) 


We impose the following compatibility conditions for the Euler boundary conditions: 


ybo{Ro) = Vp{0, Rq), ybi{Ro) = Vp{6o, Rq). (1-19) 

The boundary conditions for the remainders {u^,v^) in (1.12) - ( |1.13 ) and (1.15) are the 
no-slip conditions, and the condition in (1.17) is the stress-free condition. 


1.3 Main Result 


In order to state our main result, we must first define the norm Z, which is the Prandtl 
layer norm in which we close our nonlinear analysis: 


Definition 1.2. 


W,v\\% = 


dojdR +11 + \dR{rv)\'^r^ duidR 


.2 A 


|2 A, 


+ e^ 


u^r^+^dudR ) + 


1/9 


v^^r^+^dujdR 


1/9 


+ e^ 


uJ,fdujdR 


-he' 


7-1-1 


vf,'^r dojdR 


1/9 


( 1 . 20 ) 


where q = l + 5\ 5' arbitrarily small but positive, 7 G (0, |;). Let p be the Holder conjugate 
of q, and 0 < ^ < a < ^. Most importantly, 6 will be taken in the interval 1 — ^ < (5 < 1. 
The space Z depends on the weight 5, but we will refrain from depicting this explicitly. 


1/9 
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Theorem 1.3. Let Ub > 0 and u^{r) be a given Euler shear flow such that the derivatives 
d^u^{r),k > 1 decay exponentially. Suppose the boundary data in ( 1.12 - 1.18) are 


prescribed. Suppose that uq and ui decay exponentially fast in their arguments, that the 
compatibility conditions ( 1.19[ ) are satisfied, and that |Vbo “ ^bil ^ for small Oq. Suppose 
further that min{rtb, > 0. There exists a positive angle 9q which depends on the 

prescribed data such that for 7 G ( 0 , j) and 6 G ( 0 , 1 ) sufficiently close to 1 , the asymptotic 


expansions given in equations (1.7 - 1.9| are valid. The approximate solutions appearing in 


the expansion are those constructed in Theorems 3.1 3.5, and 3.16, and the Navier-Stokes 
remainder satisfies |< Cq. 


Corollary 1.4 (Inviscid convergence). Under the assumptions of Theorem 1.3, we have 
the following inviscid convergence: 


\U{uj,r) — Mg (r) l^r^drdio] 


+ 


<Ce'^p for 2 < p < 4, 

( 1 . 21 ) 


and for 6 arbitrarily close to 1 , and 


|t7(t(;, r) — Mg(r)|^r(irda; ) +( / / \V{uj,r)\^rdrdio ] < Ce^p for 4 < p < 00 , 

( 1 . 22 ) 


where U and V are the original Navier-Stokes flows appearing in equation (1.2). In 
we have the following convergence: 


sup \U{uj,r) 
sup |U(a;,r) 

(cj,r)Gn 


nO(r)-M°(u;,i?)|<ei+i, 
VeVpiuj,R) - yfevl{u,r)\ < e^+s. 


(1.23) 

(1.24) 


Function Space Preliminaries 

We briefly discuss the relevant function spaces in which we develop our analysis. Only basic 
definitions are given here because they are required to follow the steps of the outline below. 
The details of our functional analytic setup are presented in Section The interaction 
between the Prandtl scaling and the geometry of our domain manifests itself in the functional 
framework of our analysis for the following reason: 
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Consider an L- function, a, in Enlct coordinates. By definition, this means / J r)rdrdo < 
oo. The corresponding scaled function in Prandtl coordinates is given by u{uj,R) = u{uj,r). 

The Eulerian norm scales down to; 

11 ^ 11 l 2 (Euler) ~ J J u^rdrduj = y/l J J u^{uj^R)rdRduj ^ \fe J J u^{(jj,R)RdRduj. 


Due to the mismatch between the scaled Euler L?' norm and the actual Prandtl norm, 
we must work in a new set of function spaces (notationally depicted as || • ||* and variants 
thereof) which are natural to our problem, and build the corresponding analytic machinery 
we require to do our nonlinear analysis. Motivated by this, we define the following Prandtl- 
layer version of the norm. 


Definition 1.5. 




rdRdu. 


Applying this same analysis to the derivative operator, V = motivates the 

following definition: 

Definition 1.6. ||u||j|^i := ||w||l2 + ||V*rt ||^2 where V* = and similarly for 


where has components — 


duioj 


9 ’ 


dRR ) • Occasionally we will refer to which 


has components ^ . The corresponding weighted variants of these norms will be 

denoted with two subscripts: ||w ||^2 := // (cu, R)r^dRduj. 


Whenever we write L'p without any subscripts, this means the usual in either the 
Prandtl layer or the Euler layer, which will be clear from context. The following are the 
norms in which energy estimates will be obtained: 

Definition 1.7. ||rt||y= J J dRdu. 

J J [j'^\dR{rv)\‘^ + ev‘^r^ + dRdco. 


Definition 1.8. ||u||g := 


By the Eundamental Theorem of Calculus and Holder’s inequality: 


u^r^dRdoj < 


u^r dRdoo if u|ij=o = 0. This paired with the divergence-free condition, Uoj = 
-dfiirv), yields: J J + u^) r^dRduj < HuH^. Motivated by this, we define the following 


norm: 
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,2^-5 _l„2_^1+5 


\x — I I U r +u^r + u^r 


Definition 1.9. \W\\\ ■— 

With these definitions in hand, we detail the steps of our analysis. 


Outline of Proof 


Inserting the boundary layer expansions (1.7 - 1.9) into the scaled Navier-Stokes system 


(1.6), we obtain the following system for the Navier-Stokes remainders (for the remainder 
of this section, we replace u^,v^, by u, v, P for notational ease): 


where 


I 1 yC ve 1„ 

-UsUuj + -UsujU -h UsRV + VgUR H - VgU H - UsV + -P^ 

f e 2 3/2 „ 

- URR - Ur - kUujui + -^U - kC ' = }, 

1 ^ ry* ^ 

II 2 11 

-UsV^ + -VsujU -h VsVR + VsRV - ^UgU + -Pr 

r r r ^/e e 

\/e e 2^/e e 

— URR - Vr - kVujuj H- T^Uuj H- kV = Q, 

/y* ry*^ ry^^ ry^^ 

1 

-u^ + ve-u + vr = 0. 
r r 


f{LO,R) = —e 2 — e ')'+2 Jr VUR + 


g{io,R) = —e 2 12“ — y/eR^’^ — j -uvoj + vvr -/= — 


- ‘L 

r 


1 u‘ 
e r 


(1.25) 

(1.26) 

(1.27) 

(1.28) 
(1.29) 


Here, 12“,i?“ are the remainders from the approximate solutions Uapp,Vapp, whose precise 


definitions are given in (3.1) - (3.19). R^’P,R"’P are the linearizations of the Navier-Stokes 


remainders around the Prandtl-1 layer, which precisely are given by: 

7 ?“’^ = ^UpUuj + ^U^U P Upj^V + VpUR + ^VpU + ^UpV, 


1 


1 


„u. 


(1.30) 

(1.31) 


R^'P = -u^pV^ + -v^^u + v^VR + v^jiV - ^ul 
The NS remainders u, v satisfy the following boundary conditions: 

[itw]|t^=o = ['iiw]|R=i?o = 0) -h ruK = 0 and Pr = 2eua; on {a; = 6»o}. (1.32) 










Step I: Construction of Approximate Solutions 


We first construct the approximate solutions Uapp, Vapp such that the resulting remainder 
terms and R" are higher order in e. This involves three stages: constructing the Prandtl- 
0 layers {u^,Vp) using equations (3.1, 3.7), the Euler-1 layers using the equations 

(3.2, 3.10), and the Prandtl-1 layers (Up,Vp) using the equation (3.3). The divergence free 


conditions are enforced at each stage, and the boundary conditions are given in (1.12 


1.18). 


The method of constructing the approximate solutions is as follows: the Prandtl-0 layer 
angular velocity, Up, is constructed via a von-Mises transformation, for which the assumption 
minjub, Ug -|- > 0 is crucial. The radial velocity Vp is then obtained via the divergence 

free condition: rVp = u^. This choice creates rapid decay as 72 —)■ oo for the Prandtl-0 

layers, but as a consequence a boundary condition for Vp\R=Rp cannot be enforced. 


The second stage of the construction addresses the Euler-1 layer, (Ug, Ug) which is designed 
to correct for the normal boundary velocity of by enforcing Vp\R=Rg -f vl\R=Rg = 0. 

After passing to a vorticity formulation for Ug, this layer is obtained via standard methods 
from the second order elliptic theory. 


The last stage of the construction addresses the Prandtl-1 layer, (Up,Vp). The boundary 
conditions on {r = 72o} are UplR=Rp = —ullR=Rg and VplR=Rg = 0. These are designed such 
that UapplR=Rg = 0 and VapplR=Rg = 0. This construction relies on the positivity estimate, 
which will be discussed in Step III. 


After these three stages, we evaluate the remaining error, f f (|72“p -|- e|72^p) The 


weight of must be included because K^,R" are contained in f,g in the equations (1.28 


1.29), and accompanies f,g in the linear estimate (1.50). Interestingly, the angular 
error term arising from equation (3.4), f f |ngpr^“^, is infinite in the critical case of <5 = 1, 
which is the reason the convergence in Corollary 1.4 cannot include <5 = 1 for p < 4, which 
in turn would correspond to the usual convergence. Despite this, we make use of the 
delicate embedding theorems we prove in Section in order to recover convergence for 
p>4. 


The construction of the approximate solutions and evaluation of the resulting error 
culminates in the following: 


Theorem 1.10. Under the assumptions of Theorem 1.3 there exist approximate solutions 
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such that 


Rur^~^^ dRduj] + y/e 


\ 5 

dRduj j < e3“'^ if 0 < (^ < 1 and 


for K > 0 but arbitrarily small. Moreover, the approximate solutions satisfy the various 


estimates which appear in Theorems 3.1, 3.5, and 3.16 


Step II: Energy Estimate 


In this step, we obtain the natural energy estimate associated to the linearized system (1.25) 


(1.27) 


Theorem 1.11. 

II'^IIa ~ + <^^o|| 7’||^2 +11/11^2 +^11511^2 for (5 G [0,1] and e << 00- (1.33) 

^*,S *, 2+5 *, 2+5 


This estimate is generated by applying the multiplier to equations (1.25 


1.26). Once the weight is fixed for the angular multiplier, the divergence free condition 


(1.27) forces a loss of one factor of r. We illustrate this by multiplying the right-hand side 


of (1.25) by yielding: 


fr^+^u < / / fr^+^ + 


2^5 < 


u r 


f2r2+5 ^ 


uy 


< 

r\_/ 


J I /V+'5 + 02 J I\dR{rv)y< I J /V+'5 + 0g||u|||. (1.34) 

Thus ||u||b must appear in our estimate, which features an extra factor of r as compared 


to ||u||y 4 according to Definitions 1.7 and 1.8 The strongest weight for the radial multiplier 


which is then consistent with the presence of ||u||s is er^~^^v. 


Step III: Positivity Estimate 

In this crucial step, we estimate ||u||b in terms of ||u||a. Snch an estimate must overcome 
two difficulties. First and foremost, multiplying eqnation (|1.26) by a multiplier which is 


0(ev), as in Step II, formally results in control over e f f |VeUp, which is too weak in 
the inviscid limit. This lack of a basic order-one estimate of v is the most fundamental 
difficulty in the boundary layer theory. In the case of a rectangular geometry, Guo and 
Nguyen overcame this difficulty by using the vorticity multiplier dy ) |GN14j . 

Second, since ||u||s, which appears on the right-hand side of estimate (1.33), contains an 


extra factor of r when compared to ||u||aj the positivity estimate must recover this factor. 
This difficulty is new to our problem due to the geometry of our domain D. 
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The starting point is the following calculation, which we use in Sectionj^and throughout the 
construction of the approximate solutions. We temporarily ignore boundary contributions 
as we shall apply this calculation to functions v vanishing on relevant parts of the boundary. 

Lemma 1.12 (Positivity Calculation). 

J j r^\dR{rv)\^ < J j r^\dR{rv)\^ + J j UsUsRr + Oo\\v\\%. (1.35) 

Proof. For the sake of simplicity, we select the h = 0 case to showcase initially. The case 
for general 6 which we shall need involves controlling a few more terms, and is proved 
rigorously in Section 



jdR(rv)l =11 \dR{r—Us)\" 


2„,2 


\dR{^rul + 

Ug 


2 

„,2 ^sRi 


r V 


u 


2 ^sR 


< 


u^srII OrCA) <[ la^^)! 

\J Ro Us J J J JHg Us 


(1.37) 
2 (1.38) 


Recalling that minu^ > 0, and inserting (1.38) in (1.37) and then into (1.36) yields the 
desired estimate. 

□ 


The key calculation is that in estimate (1.38), in which we’ve used the rapid decay of Usr 
to conclude: 


sup / u‘lji{R — Ro)dR < oo 
t^G[O,0o] ^0 


(1.39) 


This will be proven rigorously in equation (4.27). Moreover, as in |GN14j . our positivity 
estimate relies on the profile tt* > 0, which in turn relies upon our assumption that Ub > 0. 
This is the reason our analysis does not treat the case of a non-rotating boundary. This 
lemma is used to prove the following; 
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Theorem 1.13 (Positivity Estimate). For 6 G [0,1], 


Ml + 


I uj=9q 


^ ^II^IIa + (^0 + Ve)||^^||s + (1 “ ^)^I|-P||r2 + " 


(1.40) 


In particular for Oq and e small enough, this establishes control of ||u||b in terms of ||um: 

,2,,. .^2noM2 , ,,.,,2 , m.M2 




/uj=eo 


euy-^ <1 


yA + ii-M\\p\\hM\\f\\h^^ye\\g\M^ ■ 

*,0 *,2+0 *,2+d 


The essential mechanism behind the positivity estimate is to capitalize on the order 1 

3r] 

free condition. We apply the multiplier MdR{—-), —ed^ji 


appearance of vr in the positive profile term in equation (|1.25|) through the divergence 


r^+^v 


Us 


)) to equations (1.25 


which is formally a weighted vorticity multiplier. The weights are designed care: 


1.26), 
ully to 

capture the ||u||b norm using the profile terms from (1.25 - 1.26). We highlight this using 
the three important profile terms below: 

j j ^uydnC-^) - j j r^u^dR{rv) + j j u9ij(ru) 


Us 


— e 


UsvMdu, 1 — 

Us 


—e 


r^\dR{rvM + ^ f j r^^dR{{rvf ), 




vu.sr^asM I r^'^BnUrvf) - f [ 


ui 


( |1.42| - | 1.44D « - 
< _ 


- 


2+5 UsRR 2 
r - V 


Summing (|1.42 - 1.44|), integrating by parts, and using Lemma (|1.12|) yields: 

r^\dR{rvM - 

r'^|( 9 ij(ru)|^ - ey J \ 


r^vl- 


(1.42) 

(1.43) 

(1.44) 

(1.45) 

(1.46) 


Once the important quantities from ||u||s have been extracted, the rest of the proof 
proceeds by estimating the remaining terms after the multiplier is applied to the equations 


(1.25) - (1.26). 


Step IV: Pressure Estimate 

The pressure term ||P||r2 must be included in Theorems 

*, 5 

effects. The choice of h = 0 for the energy estimate multiplier .n Theorem 1.11 forces the 


1.11 


and 


1.13 


due to geometric 
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pressure term to drop out whereas the choice of <5 = 1 is required by Theorem 1.13 in order 


for this term to vanish. The lack of a consistent choice of 6 which simultaneously forces 
the pressure to drop out of Theorems 


following: 


1.11 


and 


1.13 


requires us to estimate ||T’||i;, 2 ^ in the 


Theorem 1.14. For 5 G [0,1] and e << Oq, 

W^W'L'i ^ Ci{ 9 o, e)\\u\\\ + C2{Oo, e)\\v\\% + \ \f, y/eg\\\2 
*,0 *,( 

Here, Ci{9o, e) —)• 0 as either 0o —^ 0 or e —>■ 0. 


(1.47) 


We emphasize that this estimate is new in our analysis due to the presence of geometric 
effects, and therefore did not appear in [GN14| . Moreover, it is surprising that the ||u||a and 


|u||b terms appearing on the right-hand-side of Theorem 1.14 are accompanied by small 


parameters. The estimate relies on the existence of a vector field, Ai = (a(cu, R), b{u}, R)) 
such that div(Ai) ~ P, and ||Ai||j|^i < ||P||£, 2 , which is guaranteed to exist for P G by 


|Orlt98( Page 27] and the estimates we establish in Claims 6.2 - 6.5 Moreover, Ai can be 
selected to vanish on the Dirichlet portions of the boundary, {P = Rq} and {ui = 0}. Given 


this vector field, we apply the multiplier {ar^,ebr^) to (1.25 - 1.26). The weighted vector 
field is used as our multiplier in order to estimate the correct weight on the Pressure term: 


—r^a + Rnbr^ « — 
r 


+ bR)r^ ^ - j j Pdiv{Ai)r^ ^ - j j 




(1.48) 


Once ||P ||£,2 has been extracted, the rest of the proof proceeds by controlling the terms 


arising from applying the multiplier {ar^,ebr^) to (1.25 - 1.26) 


Summary of Linear Analysis in Steps II - IV: 


Putting estimates (1.33), (1.41), (1.47) together yields the full energy estimate for the 


linearized system in (1.25 - 1.27): 


W\\\ + Iklls + ll-P|l r2 + 


*’■5 Ja;= 6 »o 


(1.49) 


When paired with the divergence-free condition, we can upgrade u, from order ^/e to 
order 1: 


I Il2 I II ||2 I 110112 

mix + ll'^lls + Ik 




+ f 1 < ||/,^ 5||2 

Juj=eo 


(1.50) 
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The existence of a unique solution to the linear problem (1.25 - 1.27) is then given by an 
application of Schaefer’s fixed point theorem: 


Theorem 1.15. Let Us and Vg be the approximate solutions as defined in equations (1.7 


1.8). Then there exists a unique solution [u,u,P] to the system in (1.25 - 1.27) on the 


domain Q, together with the boundary conditions (1.32) which satisfies estimate (1.50) 
uniformly in e and small Oq. 


Step V: High Regularity Estimates 


In this step, we obtain higher regularity estimates for solutions to the problem (1.25 - 1.27). 


To do so, we rewrite the equations (1.25 - 1.27) by moving the profile-dependent terms to 
the right-hand-side: 


URR 


ye e 

Ur o ' 


^ I 

— VRR - VR - kVuju) + n 


e 

2^ 


_ ±^3/2 
^2 


+ -Pco = /, 


e 1 
Uui + -^v + -Pr = g, 


where 


/. / 1 1 "v/c "v/c 

J = J - [ -UgU^ + -UsujU + UsRV + VgUR -I- VgU ^ - UgV 

r r r r 


9 = 9 


11 2 1 

-UgVuj -f -VsulU + VgVR -b VsRV - -^UgU 

r r r yje 


(1.51) 

(1.52) 

(1.53) 

(1.54) 


From this point of view, we formally expect high regularity estimates using the standard 
theory of the Stokes equation: ||?x, u ||^2 < e~^||/, some potentially large value 

M. This is only a formal estimate, however, because the corners of fl, {ui = 0, R = Rq) and 
{u = Oq, R = Rq), obstruct the regularity of the standard Stokes problem. To account 
for this, we use the results of |Orlt98) to recover 77^/^ regularity for the solutions near the 
corners. Precisely, the main result of this section is: 

Lemma 1.16. The solutions u and v can be decomposed into u = ui + U2 and v = vi + V2, 
where U2 , V2 are supported near the corners of the domain fl in a region 712 satisfying 
{u},R) G 712 implies R — Rq < l,r — Rq < -y/e. The decomposition obeys the following 
estimates: 


+ lk2r-™,U2r™||j^3/2 


<e ^||/,^/e5||i2 


2+5 


(1.55) 


for some possibly large value of M, where the constant is independent of Oq, and for m 
arbitrarily large. 
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Step VI: Nonlinear Analysis 


In the final step, we use a contraction mapping to obtain existence and uniqueness of 


the nonlinear problem (1.25) - (1.32). Consider a sample nonlinear term, from 


equation (1.28). According to the right-hand-side of estimate (1.50), we must estimate the 

II • IIl?, 


,2+S 


norm of the nonlinearity: 




vurW 


* , 2+<5 



- 7 - 1-5 


VUR 


‘ ^ 2 +S < ^27-Hl 


y'^Pj.P-^ + ^P 


‘^‘l„q+a 


u^r 


(1.56) 


We have used Holder’s inequality and we suppose for this discussion that the technical 
parameter a satisfies ^ ^ in order to make the above inequality valid. We think of p as 

being very large and so q = is very close to 1. This calculation then motivates two 
features of our norm Z. 

First, Z must control f f and / j{^/ev)'^P for large p, together with the appropriate 
choice of weights r. Typically, the embedding yields the desired control, 

but cannot be applied in our setting for two reasons. First, the standard embeddings apply 
to integrals taken against the usual measure RdRdco and for the usual gradient operator 
V = whereas in our setting the measure is r{R)dRdoj and V is replaced by V*. 

Second, we must precisely determine the weight of r{R) that can be controlled by our 
weighted energy norms, X and B. As such, we establish the required embeddings from 
scratch. 


The second ingredient which is built into Z are high regularity quantities, because as seen 
in estimate (1.56), Z must control f f |Vw, In order to close a contraction mapping 

argument, we must in turn control these high regularity quantities (see Definition 1.2). 
The main result in this direction, which serves as the driving force behind the contraction 
mapping argument, is: 


Theorem 1.17. For u,v solutions to the system (1.25 - 1.27), there exists a > 0 such 
that if g = 1 -|- 5' and p = we have: 


e4||u^|| 2, +e'i\\vR\\j2q +e‘^\\u^ 






1^2, Te4 + i||?;^|| 2 , <||/,Ve5llL2 (1-57) 

*,0 —/3 *, 2 +d 


for all 6 such that max{l — where /3 > 0 and we can take 0 < | < a < ^. 
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The essence of the proof of Theorem 1.17 is as follows: since 2q is only slightly larger than 


2, we interpolate between the estimate in (1.50) which is uniform in e and the high 


regularity estimates in Lemma 1.16t which scale poorly in e. Again, these interpolations 


are highly sensitive to the weights r which can be controlled by ||/, y/ig\\ij 2 , and are 
also taking place in our =t=—spaces, and so must be developed from scratch. The required 
embedding theorems are proven in Section and Theorem |1.17 is proved in Section 


With Theorem 1.17 in hand, we are able to close a contraction mapping argument in the 
space Z, which we do in Section 

Theorem 1.18 (Nonlinear Existence and Uniqueness in Z). For 5 G (0, 1) sufficiently close 
to 1 there exists unique Navier-Stokes remainders {u,v,P) to the system ( |1.25 - 1.31) such 
that ||u,u ||2 < oo. 


From here, the main result in Theorem 1.3 follows immediately. 


2 Function Spaces and Embedding Theorems 

2.1 Basic Properties of ^-spaces 

Here we establish a few basic facts which will be in use throughout our analysis: 

Lemma 2.1. Holder’s Inequality for For p, g Holder conjugates, < | |u| 1^? ^ 



The third inequality above is the usual Holder inequality against the standard measure 
RdRdoj. 

□ 

Lemma 2.2. The space endowed with the norm || • is a Banach Space for 

1 < p < oo. 
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Proof. It is clear that J J uPr^ = 0 u = 0 and that J J [cuYr^ ~ J J ~ 

11 ,,V, We check the triangle inequality by ueing the cortesponding triangle inequality 

for the usual norm which corresponds to the weight RdRduj: 

u x/„,, . ,, V 






We must argue that LF^ ^{Pl) is complete under this norm. Suppose {un} is a Cauchy 
sequence Cauchy in the usual LP norm, so there exists a limit function 


u such that 


s/p 


Rl/P 


u in RP. Define u = uj^rj^, so we have: 


u, 


- uPr^dRdco = 


Un s/p 

——r 'P — u—:rr 

Rl/P Rl/P 


RdRdio —)■ 0. 


□ 


Lemma 2.3. The space is a Hilbert space, endowed with the inner product 

(n, v) = f f uvr^dRduj. 


Proof. By the Holder’s inequality (established above), the inner product is well defined 
as a mapping ^ x ^ M. Moreover, it is easy to see linearity, symmetry, and non¬ 
degeneracy of the inner product. By the previous lemma, this inner product induces a 
norm, and the space is complete with respect to this norm. □ 

In general, many properties of the usual IP will be inherited by RP ^ because the map 
T : RP g ^ RP given by T(/) = f-^^ is a linear isometry. For instance, the characterization 
of the dual space to R^ ^ follows trivially from this observation: 


Lemma 2.4. \^RP ^ j = R^ s where the superscript * denotes (as always) the dual space. 

Proof. Given a bounded linear functional / : —)• M, / o T~^ is a bounded linear 

functional IP —)■ M, and is therefore given hy f ^ J J JgRdRdoj for some g G L*?, where 

/ = Letting g = gR^/ir~'^, we readily check J J fgr^dRdu ~ J J JgRdRdoj and 


that 


\Ll , = \\9\\L<i- 


□ 
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This immediately implies reflexivity for 1 < p < oo, and thus we will be able to obtain 
weak subsequential limits from sequences bounded uniformly in spaces in the usual 
manner. We’ll need a few more facts: 


Lemma 2.5. If Un — ^ u for any 1 < p < oo and any weight k G M, a subsequence 

a.e. 

Uuk -^ U. 


Proof. Define Un = 


Un 


Rl/P 


TP 


and u = 


u 


Ri/p 


rp. Then Ur. 


u, so a subsequence u 


nk 


which immediately implies u. 


rik 


U. 


->■ U, 


□ 


Lemma 2.6 (Density of in For 1 < p < oo and any weight k G M, we have that 

is dense in 


Proof. Given an / G dehne / = 


/ 


Rl/P 


rp which is now in the usual L^. By density, 


there exists (pn f 11 \ f ~ (pnl^ RdRduj —)• 0. Now define = (pn, which 




rRp 


immediately yields (pn —^ / and moreover cpn G C“(D) because R> Rq. 


□ 


2.2 Properties of Z, I 

In this subsection, we prove the first basic property of the space Z\ 


Lemma 2.7. The space Z together with the norm ||n, defined above is a Banach space. 


Proof. Nondegeneracy and homogeneity of the || • \\z follows from the definition. The 
triangle inequality follows from applying it separately to each component, showing || • ||^ is 
a norm. We must verify completeness. We first show that weak derivatives, when they are 
elements of the space ^ for any r >2 and weight n, are unique within this space. 


Suppose we have two weak radial derivatives and in of u. Then J J {ur- u%) (pdujdR 


0 for all (p G Cf^{n). Since the support of f is compact, we also have 




u]f)(pr'^dujdR = 0 for all (p G CpP{U). Let s be the Holder conjugate to r, and select 


an arbitrary f € ^ ^)*. By Lemma 2.6 approximate / by (pn in the norm L 


' S 



— ujfifr'^dujdR 



- u r) {LI^ Jim) f nr 


K 



UR)(pnr'^ = 0 . 
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We have exchanged the L'l ^ lim and J J by using Holder’s inequality: 


- tt|) ((/)n - /) r'^dujdR 


< 


'^R 


41 


\<Pn - 


Since the right-hand side goes to zero in the above inequality, we are able to switch the 
limit and integral. Thus, has operator norm 0, and the only such element is the 0 

element, showing radial derivatives are unique within the class for any r and k. The 
choice of radial derivative as opposed to angular derivative was without loss of generality, 
so the above uniqueness result holds for angular derivative as well. 


Suppose {u"'} is Cauchy in Z. Then in particular {«"■} is Cauchy in so there exists 

Llg Llg 

limit u such that u" —u by completeness of By the same argument, v 

L2 


and w. We must verify v = Uco and w = ur. Let (p G Then: 

wpdRduj = — J J lim) = lim J J = — lim J J u^pR 

= - J J {Li s lim) = - J J upR. 


We can exchange the Ll g limit and integral again by Holder’s inequality. Since ur is 
the unique element in Ll g satisfying the above equality, we have ur = w. The identical 
argument shows = v. 


We now turn to the term. Again by completeness, there exists some limit 

j^2q 

function, w (we will abuse notation), such that w. By passing to a subsequence, 

Llg 

we can assume w. But we know from earlier that —A ur so a further 

subsequence must converge almost everywhere to ur. Since every subsequence of an a.e. 
converging sequence must also converge a.e. to the same limit function, we have w = ur. 

Since all of the weights above were done in full generality, we can repeat these arguments 
for all of the terms in the norm. This proves completeness since we have exhibited a single 
element u which serves as the Z—limit of the Cauchy sequence vL. 

□ 


Corollary 2.8. The spaces A, B, and X are individually Banach spaces. 
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2.3 Properties of Z, II: Low Regularity Embeddings 

We prove weighted embedding theorems which replace the usual ^ Sobolev em¬ 
bedding in This style of argument will be applied repeatedly in this paper. For this 
section, we suppose that u, v satisfy the boundary conditions displayed in (1.32) and satisfy 
the divergence free condition in (1.27). 

1 

Lemma 2.9. j dRduj^ < ||u||b for 2 < p < oo. 

Proof. First consider the case p = 2. In this case the exponent on the weight r is <5, and so 
we have the result by definition of ||u||b. 


Next, consider the case p = 4. We express: 


rR 


^ ^ / d^{v^)r^ / dy{v‘^r{y)^+^) 


' Ro 


r-K 


rR 


10 


vVuiT^ / + / vvujr^ / v‘^^/er{y)^ 


' Ro 


I Ro 


rdo roo 

< / \vv^\r^ 

Jo J 


poo pOq poo 

/ + y/e / \vVi^\r^ / 

JRa Jo Jr 


( 2 . 1 ) 


Integrating both sides in dujdR and applying Holder yields: 

1/2 


4 l-h25 < 


V r 




,.2„<5 


Ro 


\ 1/2 


+ \/e 




,.2„<5 


( 2 . 2 ) 


We now have to distribute among the right hand side of the inequality, which 

yields the desired result. For p G (2,4), we interpolate: 


„ p 1 I 5p \ P 

vPr2-^+^dRdw\ = 


|ur2 + 2|Pr ^dRduj 

1 I <5 Q -j ^ 

|ur 2 ^ 2 1 r~ dRdu; j 


< 

< \\v\\b- 

Now for p > 4 we can proceed inductively via the calculation: 


1 I S A -t ^ 

r~ dRdcj 


1-0 

4 


(2.3) 


£_1 I E E_i4_l£ E E4_1 e 

yPjr.2 ^^2 = i; 27’4 ^■^4t!2r4'^4 = 


E_i E_i_|_1e / E_i / nE_|_EE 

V2 Vu)r4, ^ V2 vjir[y)* * 

^0 J Ro 


rR 


, E_l^£ 
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2-1 , ^2-1+«£ 




+ Ve y VujT ^.' 


p / nE-|_^£_i 


'Ro 


v2r{y)‘ 


(2.4) 


Taking absolute value, integrating, and using Holder yields: 

1/2 / r r \ 1/2 


E_i+«E| . 

y^'f 2 ^ 2 <C 




vlr^ 




+ ve 


^p-2^2^-l+f(p-2) 


^2 5 
*^0;' 


p Pt§p_-i 

y 2 y 4~'~ 4 


(2.5) 


If p is an even integer, the absolute values on the left-hand side can be removed. We 
therefore establish the inequality for even integers successively starting at p = 6 (since p = 4 
has been computed directly), and then interpolate in between. 

□ 


Sp 1 \ 2 

yPy 2 r 2 < ||u||x for 4 < p < oo and 


1 

Sp \ P 


u^r 2 ^ ||^i||x foi* 


Lemma 2.10. 

2 <p < 4. 

Proof. For p = 2, we have f f < ||rt||x t>y dehnition of the norm. We compute the 
case p = 4: 

rR 1 i-R 

uuRr{y)2+^ + Ve 

I Ro J Ro 


«^r2'5+2 = 


< 


[ uuu,r^ ( [ uuRr{y)2+^ + ^/e [ u^r{yf 2^ 

0 \J i?o Rq / 

pOq poo ^ p9o poo ^ 

/ / \uuR\r2^^ + / \uuui\r^ f u^r^~2. (2.6) 

Jo Jro Jo 


' Rq 0 J Rq 

Integrating both sides over dudR and applying Holder’s inequality yields: 




u r 


+ 






1/2 


1/2 


.2„& 


.2 ...<5 


1/2 


1/2 




1/2 




1/2 


2 


^ l|ii|lx- 

We now interpolate for p G (2,4): 


(2.7) 


1 e 

J j |rtr 2 ~ I J 


u‘^r‘^^dRdu?\ 


1-6 

4 


< l|ii||x- 
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Once the above estimate for p G (2,4) has been established, the desired estimate can be 
inductively established via: 


yPr2 


Sp 


p ^ p Sp I 1 

= u2r i w^r -i 



p_i ^ 

u2 Ui^r i 



p Sp I 1 

dR{u2r{y)^^2) 


E_1 §P 
M2 r 4 


f-R 


2-1 

M2 


' Ro 


URr{y) 4' 


1 ^ f‘^ E_1 §£. 

2 + Ve / 4 

Jo 



1 

2 . 


( 2 . 8 ) 


Taking absolute values and applying Holder’s inequality yields: 

I jyF-\up-^) {^j juyy" y 
+v~el Juyyyi juP-yip-y"^". ( 2 . 9 ) 

For p > 4 all of the quantities in the above estimate are inductively controlled by powers 
of ||m||x- 

□ 



Remark 2.11. This argument is reminiscent of the proof of the classical Gagliardo- 
Nirenberg-Sobolev inequality. This method can be used to yield a direct proof of the 
standard ^ LP embedding in by replacing the weights r by R. The advantages of 
the direct approach above is the avoidance of defining the fractional Sobolev spaces 
and consequently the avoidance of appealing to the Fourier Transform. Indeed, in the 
classical case, one must argue ^ “—)• LP for 0 < s < 1 and 2 < p < oo because the 

Sobolev exponent is critical. The drawbacks are that this method must take into account 
the behavior of u on the boundary dLl, and that it doesn’t directly apply to more complex 
domains. 

We now prove an embedding of the type Z L °°, from which Corollary |1.4| follows 
directly: 

Lemma 2.12. Given d' > 0, let g = 1 + 6'. For < 5 < 1, 

£2 + 49 I _|_ g2+4<j + 2 < I u| 1^. (2.10) 
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Proof. Since 2q = 2(1 + S') > 2, by Money’s Inequality: 

\u{Co,R)\ < \\Vu\\l2<i + \ \u\\l2, = (^j j l^ul'^mdRdujy" + (^j j \u\‘^'^RdRdi 


\UJ 


2q 


( 2 . 11 ) 


Multiplying by e*<i yields: 


g49 + 2 \u{ljJ, i2)| <62 


\Vu\^'^rdRdLo\ + 


j_ 

2q 


\u\^'^rdRduj 


2q 


< 62 


+62 


\u\‘^‘irdRduj] 


2q 


< 


\‘U\\z- 

( 2 . 12 ) 


where we have used |Vn| = | I ~ I I = |V*ri|. The condition < <5 


ensures ||w|L 2 g < ||tt||r 2 q < \\u\\z by Lemma 


,5q 


2.10 


The proof for v works identically, 


where the extra factor of 6 2 is required as the \\v\\z contains 6 2<2 J v^r 2 p^ 


□ 


Remark 2.13. Note that the condition in Dehnition 
condition < <5 < 1. 


1.2 


1 — ^ < (5 < 1, implies the 


2.4 Properties of Z, III: High Regularity Embeddings 

In this subsection, we provide careful estimates which will yield control of the high regularity 
quantities appearing in || • \ \z. Throughout this section, u,v are assumed to satisfy the 


boundary conditions displayed in (1.32) and the divergence free condition in (1.27). 


Lemma 2.14. Let <5 G [|, 1]. There exists a 5' > 0 such that if g = 1 + (5^, then 


IMI 


r 4 

*,2+^ 


1 

I 2 I 

Ixl 


I zt2 
-"*,2+5 


; and ||uij||^4 


Lf ll‘'l+2 

-^*,2+5 


(2.13) 


for all 0 < a < Moreover, a can be selected such that - = htT < q- < ^ where p is 
the Holder conjugate of q by taking S' small enough. 

Proof. After noticing that u = = 0 on the boundary {w = 0} and vr = ^{ucj + y/ev) = 0 

on the boundary {R = i?o}) the u and v estimates follow in an identical manner, so we 
focus on u. We express: 

= u%r‘^^°'' = u|r5+“'n|ri. (2.14) 
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Since S > ^, the quantity uj^r^ is integrable and lies in the Sobolev space for 

each hxed u) by the definition of H'^ 2 + 5 ^ implying that this quantity decays at 00. This 
enables us to write: 

POO POO POO 

u|r^/2 = _ / dR{u%{y)r^^‘^)dy< \uRURR\r^^^dy + y/e \uR\r^^'^dy. (2.15) 

Jr J Rq J Rq 

We also note that u = rt/j = 0 on the boundary a; = 0. Thus, we are able to write: 


,2 ^k+a' _ 


URURujV 


1 +"' < 


\URURj_,\r 


1/2+a' 


(2.16) 


Multiplying the previous two inequalities, integrating and applying Holder’s inequality 
yields: 




7/2 rl+“' 
UrT 


2 2 
'’J-rr''' 


,,2 a' 

URu^r 


II Il2 II ||2 

where we use that = 2^ < 5. Taking fourth roots yields the result. 


(2.17) 


1 1 

Lemma 2.15. There exists a 5' such that for q = 1 + 5\ ||wa;||L 4 < ||rt||xll^ll ?r 2 

^*,2+S 


S G 1 ]. 


Proof. Writing = u^r and recalling that u = Uuj = 0 on R = Rq, enables 

us to write: 


u^.,r^ = 


pR pR pR 

/ dR{ulr^)dy = / r^u^u^R + ^/e ulr^~^. (2.18) 

J Rci J Rci J Rc\ 


Using the divergence free condition, u^o = —dR{rv), we have Uoj = 0 on cj = 0. Therefore 
we can write: 


u^r = / UujUujujV = / Uuir^UuiuiT 


§_ _35 

’ 2 ?/„,r 2 


Multiplying the two equalities above together, taking absolute values, and applying Holder 
yields: 


.2 „-3(5 


where we have used 5 < <5 < 1 —3/5 < —2 + 5. 


uIrT^] <\\u\\\\\u\\%2 (2.19) 
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Lemma 2.16. There exists 5' such that for q = 1 + 5', we have -v/e||'(;[j||£^4 < 

II'^^IIbII^II ^2 /5 > 0 and 1 — ^ < 5 <1. 

^*,2+5 ^ 

Proof. Temporarily writing /3' = ^, we proceed to write: v‘^r~^' = v‘^r~^'r~^v‘^r^. Using 
that u = = 0 on the boundary R = Rq, we can write 


= 


rR 


rR 


rR 


' Ro 


dy{vlr^)dy = / VujVujRV^ + y/e r 


I Ro 


I Ro 


Next, we recall the boundary conditions at a; = 0o are v^j = —ur —)■ v^{9o,R) = 

J.2 ^ 

-wUr{9o, R)"^■ As such, we write: 

„ 2 —/ 3'-(5 r9o 

vlr~^ r~^ = -^— ur{9o, Rf + / v^Vuiujr~^~^ . 


Taking absolute values and multiplying the previous two equalities together yields: 

^ ^ 2 -/ 3'-(5 roo rdo ^ roo 

vtr~^ < -s— UR{9o,Rf / \v^VujR\r^dy + / \vujVuiuj\r~^~'^ du / \vojV^R\r^dy 




' Ro 


+ 


poo 

UR{9Q,Rfyfe / 
jRo 


v^r ^+'^+ / \v^Vujuj\r 

Jo 


10 

rOo 


' Rq 
roo 

/ vlr-^+^dy. 


I Ro 


Integrating over duj and dR yields: 




1/2 


+ 


+ 




1/2 


Vuuir 


-35-2/3' 


1/2 


vlr^ 


1/2 


r3/2 


roo 

Jro 




+ \/e 


,,2„5 


2-1+5 

^(Jj‘ 


2 -35-2/3' 


^2-1+5 


,,2 „5 

^ojR^ 


1/2 


( 2 . 20 ) 


As u = Ur on the boundary 0 = 0, we can write ur{9o, R)= ur^ +■ ur{9q, RY 

_ Jo 

2 

/ uYuj- Inserting this into (2.20): 


< 

rs_/ 


1 poo p9o 

^ ^ / [ r 

^ Rq J 0 


2-/3'-<5„2 


i?a; 


v^,r^ 


1/2 


,,2 „5 

^uiRT 


1/2 
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+ 


vlr^ 


1/2 


Vux.r 


- 35 - 2 / 3 ' 


1/2 


vlr^ 


1/2 


7,2 ^<5 


1/2 


+ 


1 

+ V~e 


2-0'-S 2 
^ ^Ru 






2 - 35 - 2 / 3 ' 

^(Jjio' 




( 2 . 21 ) 


We multiply estimate (2.21) by e^, and according to the weights we are able to estimate 
for the H^ 2+5 terms, we require: 


2-(3' -6 <6 ^ 2-/3'<25 


and 


-35 -2(3' <5-2 2-213' < 45 

The lemma has been proved. 


2 “ 


1 W<,, 
2 2 - 


( 2 . 22 ) 

( 2 . 23 ) 

□ 


Remark 2.17. The above estimate for is the most delicate of the high-order estimates 
as it relies on the stress-free boundary condition placed at the boundary {oj = 6 * 0 }. For 
our purposes in Section we will take /3 = ^ in which case the valid interval for 5 is 


3 Construction of Approximate Solutions 

In this section, we construct the approximate solutions in the expansion ([^[g, and 
estimate the corresponding errors and R". First, we record the errors R^,K", which are 


obtained by inserting the expansion (1.7 - 1.9) into the scaled NS system (1.6): 
Angular Error, 

e° order error: ^(u° -h u°)u°^ -h (u° 


e^/^order error, Euler: + + 

r 


u2ul 


^''^order error, BL: -{ul -h Mp)n° -h -h —{ul^ + ulj + ul^) 


1 7^^ 1 

+ + V^U^rVp + VpUpH H ^ -|- -Up{Vp + Ug) H ^ ^‘pRR “ ~''^pR' 


( 3 . 1 ) 

( 3 . 2 ) 

( 3 . 3 ) 
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order error; -{u], + u^)duj{u], + Up) + (r;® + Vg)Ugj. + VpUpj^ 

^p)(^p + ^e) + p^Puj ~ '^err ~ ^(^er + ’^pp) 

order error: + -{ul + u^vl - - -ul^ - -(^Lo; + i^Lj 


(3.4) 


1 2 

Z9 (^e 4~ ^p) Z? ^ecj)) 


order error; -^r 


.1 


r 


2 '^pcp- 


(3.5) 


(3.6) 


After splitting into Euler and Boundary Layer variables, we have the following errors for 
the Radial component; 


Radial Errors, K": 

e~^ order error, BL: Ppji; 

order error, Euler: — 


(4?. 


e order error, BL: Ppp — 


«) 


r 

0^2 

p, 


_ 2 _A_p. 


n , 1 ril 1 

e order error, Euler: P^^ -\ - — 2 ——; 


e° order error, BL: Ppp + - {u^v^ + (i;° + Ug)) + {Vp + vl)Vpp 

~ ~ ((^e + Up)'^p + ^p^e) “ ^pRR'i 

order error, Euler: -u\v\^^ P ulvlr -('*^e)^! 

order error, BL: ^ (n^u°^ + Mp5a;(u° + u^)) + ^((u° + u°)u^^) + v^v^p 

+ 'Vpi^lr + vIr) + vKvIp) - ^{^ulul + (u^f) - ^v°p + - vIpp, 

11 2 

order error, Euler: - - -vl^ - + -|; 


order error, BL: -ulvl^ + -ulv^^ + vlvl^ + v^v^p - -vlp - + ^ul 


■^p^puj 

1 

+ -^v, 


Wp l/g,. 


^2 


2 


0 . 

2 "P’ 


order error:- ^vl, 


2 pujcj' 


(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 
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3.1 Prandtl-0 Layer 


We obtain the Prandtl-0 layer equations from the e® order angular error, equation (3.1 1 , and 
the order radial error, equation (3.7) . We also enforce the divergence free condition. 


Thus, after dropping the subscripts, the Prandtl-0 layer equations are the following: 

(Ug -h u)uui + r{vl -h v)uR + Puj = ruRR, (3.20) 

Uco + y/ev + rvR = 0 , 

Pr = 0 . 

The boundary conditions we take are: u{u},Ro) = Ub—Ue] u{0,R) = uo{R); and v{uj,Ro) = 


—vl{uj,Ro)- As will be shown rigorously in Theorem 3.1 the Prandtl-0 profiles u,v decay 
sufficiently rapidly, and so evaluating the equation above at i? = oo yields P^j = 0. This 
implies the Prandtl pressure is constant when coupled with Pr = 0. We rewrite the first 
equation as: 

{Ue + u)uui+ (rv + Rovliu, Rq)) Ur = RqUrr + 60 + 61 + 62 (3.21) 

=> (Ug + u)u^ + {rv - Rov{uj, Rq)) ur = RqUrr + 60 + 61 + 62. 


Here we have defined: 

60 := \/€{R - Ro)u°r{'i")uLo + yA{R - Ro)dr{rvl{uj,r))uR, 

61 := \/e{R — Ro)urr, 


(3.22) 

(3.23) 


rR rR 

62 :=u^^/e {u^^^{r{r])) - u^^j.{r)) dr] + UR^/e dr{r{9)vl{u!,r{9))) - drirvl{u},r))de 

J Rq J Rq 

pR pr] pR pr] 

= / / Ug^j.{r{9))d9dr] + euR / / d‘^{r{9)v\{uj,r{9)))d9dr]. (3.24) 

J Rq Jr J Rq J r 

62 is high order in e, as will be demonstrated in a later section, eo and ei will be put into the 
Prandtl-1 layer, and so we are left with solving [ue + u)uu + {rv — Rqv{uj, Ro))ur = Rqurr 
together with the divergence free condition and the boundary conditions described above. 

poo poo 

To satisfy the divergence free condition, we take rv{uj,R) = — dR{rv) = / Ucj, which 

Jr Jr 

ensures the profile v decays at 00 . We have the following: 

Theorem 3.1. Suppose mm{ub,Ue + uo} > cq > 0. Then for 9o sufficiently small, there 
exists a unique solution Up{u},R) such that: 

sup -b \\R'^^‘^d!j,dRU^\\L2]o,eo),L+R+) < C. (3.25) 

[O,0o] 
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Corollary 3.2. 

sup < C. (3.26) 

[O,0o] 

Proof. The proof follows exactly as in [GN14] . using the von-Mises transformation. Let 

fR 

Ue := We(-^o), and define r]{oj, R)= {ue + Up{LO, y))dy, a(w, rj) = Ue + Up{u!, R{r])), where 

jRo 

we’ve used that for each fixed w, the transformation {uj,R) —)• {co,r]{co, R)) is invertible by 
appealing to the maximum principle. We now compute: 

riR = Ue + Up{u},R) = a, 

rR f‘R 

Vi.o= Upuj = - dR{r{y)v)dy = Rov{uj,Ro) - rv{uj,R), 

J Ro J Rq 

+ = auj + ap{Rov{uj, Rq) - rv(uj, R)), 

UpR — OipTjR — OiOp. 


Inserting these identities into our equation; 


aOij + aap{RQv{uj, Rq) — rv(uj, R)) + (rv — Rov{cj, Ro))aap = Ro [aa^ + a^Upp) (3.27) 

— Rq {o.Olp')^ . 


On the parabolic boundary of our domain, a{0,r]) = Ue + Up{0, R{r])), and 0 ( 0 ;, 0) = 
Ue + np(a;,0) = Ub, both of which are strictly positive by assumption. Using the Parabolic 
maximum principle, a > C > 0 for some constant C, and therefore our equation is 


nondegenerate. The equation (3.27) along with the boundary conditions is identical (apart 
from the constant Ro) to that in [GN14) . and so the rest of the proof follows in the same 
manner. 

□ 


Co, Cl will be solved for in the Prandtl-1 layer, and so the contribution to the angular error 


is 62 , given in (3.24). 


3.2 Euler-1 Layer 


The equations for the Euler-1 Layer arise from equations (3.2) and (3.10) together with 
the divergence free condition. We drop the subscript for u, v and P within this section. 
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with the understanding that the unknowns appearing are that of the Euler-1 layer. The 
equations read: 


^21 1 

—Vui - v2u -I- Pr = 0, —u -I— -V + -Puj = 0, u^ + v + rvr = 0. (3.28) 


As described in ( 1 . 12 ) - (1.18), the boundary conditions are as follows: 


v{uj,Ro) =-Vp{uj,Ro), v{0,r) = Vboir), v{9o,r) = Vbi{r). 


(3.29) 


We go to the vorticity formulation in order to eliminate the pressure term: 

■0 2 


0 = dr + ru^v + U^V + - dujC^Vuj - + Pr) 


= UprUui + UpUujr + + TUp^rV + TUlrVr + UlVr + U^rV + 


Vujuj T UgUui Pr( 


,0 


M" 2 

= u'^Uujr + 'Tvy^rr'^ + - 

vP 2 

= Ug {- 2 Vr - rVrr) + 'TvPerr'^ + -+ -U^U^j 


= —UgrAu -|- {ru^r + ^er) U + U^ ^ —2u, 


— 2 vr -u ) := ttgPu, 


(3.30) 


where we have defined the linear operator L through equation (3.30) for ease of notation, 
and since ttg > 0, 0 = Lu if and only if 0 = u^Lv. Define the following boundary layer 
corrector: 

uj \ u°(w, 0 )^ 


B{uj,r) = ( 1 - W ) "^'7"V bo(r) + ^ ' ^l vbi{r). 


9oJ uO(0,0) ' 0o^°(0o,O)'“^^'^' 

Due to the compatibility conditions, B satisfies the same boundary conditions as v. Define 


P(w, r) = —rAB + 


ru, 


) P + ( -2Br --B]=LB 

Ug Ug 


(3.32) 


Since B and all of its derivatives decay exponentially fast, and since by assumption 
\dr{Vbo — Vbi)| < 6*0 we have \ \{r)^F\\^k,p < C where C independent of 9q. Next, for y a 
cutoff function supported on [ 0 , 1 ], define 


p6(w, r) = x(^)P(0, r) + x( ‘^° ^ ^ )F{9o, r), for e << 6 *o. 


(3.33) 
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Since each differentiation of the cutoff function gives e it is easy to see that 11 {r)'^d^Ei,\\Li < 
Ce We also record for future use that = 0- Consider w = v — B, then 


Lw = Lv — LB. In (3.67), v solves Lv = Et instead of Lu = 0 and the error made by this 


is accounted for in (3.72). Therefore 


Lw = —rAw+ 1 r 


n + 1 ^ + ( - w ] = El, — E := f; w = 0 on dCl (3.34) 

wi Wi J V r ' 


Since B is arbitrarily high regularity, obtaining estimates for w suffices to obtain estimates 
for V. Eb = F on {oj = 0, Oq} implies f\uj=o,eo = 0- 

3.2.1 Estimates 


Multiplying (3.34) by w and integrating by parts yields: 


-r{Wrr - — - -^)w - 


rwl + '^< N{5) 


2wWr = 


fw H— w'^ — 


f r + (5-hi \u^rr + 


r 

0 „m|| 

er' M oo 


<N{5) / / f^T + {5 + et,) 


w: 


ru 


err , \ 2 

u2 u2' 


w 

r 


(3.35) 


We have used the rapid decay of the derivatives of Ug. For 9q sufficiently small, we obtain 


w, 


rwi + ^+ —< 


w 


11 ''^here the constant does not depend on small 9q. 


r J J r 

We obtain weighted estimates, for n > 1 by testing the above equation against 

r'^w: 


^w‘l,= I I fwr'^ — j / + ^)r"'t(; + 2 j f WrWr"^+ 2 / / w~r 


2„n-l 


< 

r\_/ 


f2^n+l + 5 / / + eo 




w. 


2^n+l ^ yj2^n-l ^ < / / f2^n+l 


(3.36) 


With estimates in hand, we can establish existence and uniqueness. Again to ease 
notation, we define 


,0 


,0 


Lw ■.= -rAw = f- + + (^^1 +‘^^r ■=g- 


(3.37) 
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Lemma 3.3. There exists a unique solution w G to Lw = / in n, rr'Ian = 0- 

Proof. First, consider the following problem posed on the bounded domain Oat = {w G 
(0, 9q), R G (-Rq) Rq + R^)}- 

= 5 on = 0> (3-38) 

where OQn includes an additional boundary component, R = {R = Rq + A^}. It is clear 
that obeys estimates given above uniformly in N, so once each has been 
constructed, we can send N —)• oo. We hrst note the following version of the fundamental 
positivity estimate: 



Define the bilinear forms K[w,ip] := f fVw ■ Viprdrdoj + f f 
K[w,ip] := K[w,(p] — f f Using the positivity estimate, it is clear that K satishes 

the hypothesis of Lax-Milgram, but due to lack of coercivity we cannot directly apply 
Lax-Milgram to K. 


Note K[w,ip\ = f f fp K[w,(p] = f f fp + f f w = T~^{f + ^w) where 

T~^ is the solution operator to K w — T~^ is compact and 

self-adjoint, so the Fredholm alternative applied to the operator / — T~^{2-) enables us to 
conclude there either exists a unique solution to the original problem K[w,(p] = f f fip or 
there exists a nontrivial kernel. The latter option is ruled out by the estimates given 
above. □ 


3.2.2 Estimates 


Our starting point is equation (3.37). The boundary layer corrector is defined such that: 

Wujuj = 0 on {w = 0, 0o}- (3-41) 
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Differentiating (3.37) in io gives the third-order eqnation with Neumann boundary condi¬ 
tions: 

T^Wijj = ~ 0- (3.42) 

Applying the multiplier and noting that g\u)=ofio = 0 yields: 

,,2 




QujWijj — 


gWujLo < \\rw\\Hi + \\rf\\L2 + s 


w: 


(3.43) 


Applying the weighted multiplier r^Wu, gives weighted estimates inductively. By using 


(3.37) to express Wrr in terms of the rest, we have the full estimate: 

\\r'^w\\H2 < C. 


(3.44) 


We differentiate the equation (3.42) again in w, giving the Dirichlet problem for 

- rAwujuj = gojLj, w^ujlao, = 0. (3.45) 


Multiplying (3.45) by Wujuj gives: 


2 , 
rw^LOT + 


< Ik 


|2 


^ C ll/cjr 11^2 


(3.46) 


Weighted estimates are obtained inductively via the multiplier The estimate for 


Wrruj is obtained via equation (3.42), and the estimate for Wrrr is obtained by differentiating 


equation (|3.37|) in r to write Wrrr in terms of the other third order terms. This gives: 

(3.47) 


Kw\\H3<e-^/\ 


Evaluating (3.45) at cn = O,0o gives the boundary condition 

= on {uj = 0 , 00 }. 


Wa 


Differentiating (3.45) gives the fifth-order equation: 

V - gLOLOUJl 

to which we apply the multiplier 


—rw. 


rrujULo'^LJUu} 




'(jJUJLJLOU) '^LJLOUJ 


TW^ -\- 


wt. 


(3.48) 


(3.49) 


(3.50) 
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For the boundary terms, we estimate: 



1 

“ Wcjijjuj 'iUijcjttia; 

r 


uj=6o 
LO = 0 


(3.51) 




u}=eo 

LJ=0 


- / '^UlUJLjd^UJUjdT 


^ <5 11 ^ Wutuiuiu] 11 x,2 


u}=eo 

LJ = 0 

+ I I If I I t73 + C* 


w, 


'UJtjJLOLoF'ujU) 


On the right-hand side of (3.49), we have up to harmless factors: 


< - [ [ 


oj=9o 

uj=0 


(3.52) 

(3.53) 

(3.54) 


T=,,l 


F 511'JUojojcjr I \i,2 F I Ire I \ff3 F (511 rUijijijj. I \j^2 4“ ^11 ^ \ \j^2 F € ■ (3.55) 

Again, these calculations may be repeated using the weighted multiplier, to 

obtain weighted estimates. Putting the above calculations together gives: 


F r"+^re^ < e“^ 

' ^iJCJUJUJ ' ' ^rCJLJUJ rsj ^ 


(3.56) 


The estimate for Wrru)u) can be obtained from (3.45), the estimate for Wrrru) may be 


obtained by differentiating (3.42) in r, and finally the estimate for Wrrrr may be obtained 


by differentiating (3.37) twice in r, ultimately yielding: 

< .-3/2 


(3.57) 


3.2.3 Estimates 


estimates are obtained using the framework of Agmon-Douglis-Nirenberg, |ADN59) . 
where fe < 4 and q G (l,oo). To do so, cover the interior of the boundary {r = Rq} 
using one open set, W, the boundaries {a; = 0 , 0 o} using and the interior of the 

domain using LC^. Let denote the partition of unity associated to this covering, and 

yja,b,c,d ._ ^a,b,c,dy^_ yja,b,c,d equation: 


— rAw^ = —rw'^j. — — 








X 


— 2 rijj^Wr — - — 

r 


= VfF 

2 


^err , ^er 

0 uO 


u'i 


w^F 


2w^ F -w 
r 


X 


=: /^, 


X = a, b, c, d. (3.58) 
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The estimates for w'^ follow from the standard interior W^’'^ estimates: 

~ 11/11^9(0) + l|w^llvKl.9(0) < II/IIl 9(0) + ||'fi’||H2(Q) < ^(^o), 


^ \\fu}\\Li{n) + ||^^||iy2,9(o) < C{6o) 


-i+i 

- <3 

5 

--2+7 


(3.59) 

(3.60) 

(3.61) 


C{6o) is a constant that could depend poorly on 6 * 0 . The weighted estimates, 
are obtained by using the weighted estimates, is supported away from the corners 
of the domain, and so we can repeat a similar analysis as for remaining cognizant 
of the boundary condition d^w^\r=Rf^ = 0 for A: > 0. It remains to estimate w°‘'^ which 
follow by taking odd angular extensions across the boundaries {w = 0 } and {a; = 0 o}- For 
concreteness, we proceed to treat the w°' case, with the ufi estimate being identical. Define: 

w°‘{u}, r) = r) for uj G {—Oq, 0), w°‘ = w°‘ for to G (0, Oq) 

Applying a cutoff function: 


= x(‘^)^^“, supp(x) c (-y, y) X (i?o,Oo) 

ensures that w‘^ satisfies the boundary-value problem: 

-rA^“ = -rA (x(a;)«i“) = x(^)r + — 

r r 


(3.62) 

(3.63) 

w°-\r=Ro = w°'L=-eo,eo = 0, (3-64) 

where is the odd angular extension of f^. Moreover, w°' G W^’'^ whenever w'^ G 
by the condition w = = 0 at w = 0. Applying the standard estimates to the 


boundary-value problem in (3.63) gives: 

||'4)“||w2.9 < C{6q). 


(3.65) 


We can differentiate the equation (3.63) in u twice as vanishes on a neighborhood of 


{a; = —0o,Oo}, and repeatedly apply the estimates. The full estimate is then 


recovered using the same procedure as in estimate (3.47), and the full estimate on 


is recovered using the same procedure as in estimate (3.57). Combined with the estimates 
on we have established: 


Lemma 3.4 (IT^’'*' estimates, q G (l,oo)). 


I^IIw2.9 ^ C{6o), 


+ 1 + 3,9 < < 17 ( 6 * 0)6 ^+ 9 , 


+ 1+4.9 ^ <17(00)6 


-2+1 


(3.66) 


where C{6o) depends poorly on Oq. 
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3.2.4 Construction of Euler-1 Layers 

The Euler-1 layers are defined to solve: 


- u^rAvl + {ru%^ + vl + (^-2vr - ‘^vl^ = (3.67) 

(3.68) 


ul{uj,r) = ul{0,r) - / dr{rvl)de. 
Jo 


The pressure is defined to solve equation with in (3.28 ) exactly: Pg^(a;,r) = 

— The error made in the P^j equation in (3.28) is estimated as: 


+ rv2.^v -t- vP„v + Pu) 


(3.69) 


.oo 2 
-1 
r 


= + rug^u -|- ttgU — / -u^Uui + u^Eh + u^{rvrr) + 3ttgU,. -|- ( —^ — u^j. | v. 


2^ 

r 


By direct computation: 

dr{UgUuj) = —Ug^V — rU^^Vr — 2u^Vr — rU^Vrr, (3.70) 

dr{rUgj.v) = + ™err''^ + rU^^Vr, 

dr{UgV) = Ug^U -|- U^Vr- 

Each of the three terms above is known to be in and therefore decay at infinity. We 


can write: 


UgUoj + ru%v + UgV = — 


dr{u^Uuj + ru^^v + u^v). 


(3.71) 


Matching these terms with those in the integral in equation (3.69), the only term remaining 


is: 


u^{6)Eh{uj, 6)d6. 


(3.72) 


This represents the second contribution to the angular error. The estimates for the Euler-1 
layer are summarized: 


Theorem 3.5 (Euler-1 Profile Estimates). 


ir-^Ugqioo + < C, ||rX|lH3 < Ce-2, ||rX|lH4 < Ce-2- (3.73) 

IKvlWw^,, < C(0o), \\r"vl\\w3,, < C{6o)e-^+'^, ||r-’^Ugi||w4,. < C(0o)e""+% (3.74) 

where C{9q) could depend poorly on Oq, for q = (l,oo). By definition of u^, we have: 

lk"''“ellHl < C, llr^Mglloo < C. (3.75) 
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Proof. Only the uniform bound on Ug must be proven. To do so, we use the following: 


\ul{co,r)f = 


dr{rv){9, r)d6 


< [ \drirv)\‘^ {e,r)de 
Jo 


9 n 

sup |Mg(a;,r)| < / \dr{rv)\ (0,r)d0 => ||Ug||^ < sup / \dr{rv)\‘^{9,r)d9. (3.76) 

[0,6o] Jo r£[Ro,oo)Jo 

Calling ip{r) = / \dr{rv)f d9, the Sobolev embedding in gives: 


sup W{r)\<\W\\L^ + \\drip\\L^ < f j \dr{rv)\^ + [ j \drr{rv)\^ ^Wr'^vlWn^ <C. 
r£[RQ,oo) J J J J 

(3.77) 


We can proceed inductively to obtain weighted estimates on ||r”Ug||oo- These estimates 
are independent of small 9q. 

□ 


3.3 Prandtl-1 Layer 

3.3.1 Galerkin Formulation and a-Priori Estimates 

In this subsection, we solve for the Prandtl-1 layers, Up,Vp. For this subsection, we drop 


the subscripts on Up,Vp. The starting point is a modification of equation (3.3): 


(Ug -I- u)up^ + (Ug + Up)u^ + Upul^ + rUpUgg -h r{vp + vI)ur + rUpRV+ (3.78) 
(Ug + Up)Vp + UpVl + + Ef) + El + y/eu^j.v = turr + U^r, 

where we have included Eq and Pi, the contributions from the Prandtl-0 layer construction: 

Eq-.= {R-Rii)ul.,{r)ul^ + {R-RQ)dr{rvl{u:,r))ulR, Ei ■.= {R - Ro)uIrr. (3.79) 

The boundary conditions are: Up(a;,Po) = 0]Up{uj,Ro) = —ul{LO,Ro)]Up{0,R) = ui{R). In 
order to solve the order error in the radial equation, we take: 


/•O 

Ph = 

Jr 


~ «)2 


n 

, , -I- 2v2—^dt 

rit) V(t) 


pi _ 
^Plj — 


vRvR, 


2^1^ (3.80) 

Ir r{t) V(t) ^ ^ 


By the rapid decay of the Up terms, we have |d“Pp^^| < R "for arbitrarily large n and for 

OO., TJ\ _ I n,0( 


any multi-index a. Let u^(lu, R) = ttg(r) -|- Up{w, R). We define 

r — U^Up^ UpU^co ^ + ^R W’o;) 


(3.81) 
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which gives: 


u^u + u^Uoj + r{vp + vl)uii + ru%v - turr = F - Ei - Eq. ( 3 . 82 ) 

We take Or of the above equation, use the divergence-free condition, and divide by yP to 
obtain: 

-dRR{rv) -b -^ruRRV - ■^dR{ruRR) = ^ {Fr - {Eq + Ei)r) ( 3 . 83 ) 

{u^RU -b u^UR + Or (r(up -b vI)ur)) := G. 

Now we take of the above equation and again use the divergence free condition: 

-Orr (rvuj) + ^ruRRVuj + ^Or {rdRRR{rv)) v = - du,{^rvPRR)v ( 3 . 84 ) 

+du> {ruRR) ■ 


With an eye towards obtaining a weak formulation of the equation together with a-priori 
estimates, we rewrite the third term on the left-hand side as follows: 


\dR [rd%{rv)) =rd\{^d'^R{rv)) + ^e^d%{rv) + rd\{^)d\{rv) 

LL LL LL LL 

- 2rdR{^)d%{rv). 


(3.85) 


Therefore, our equation now becomes: 


- dRRirvuj) + ^ru%RVu; + rd]i{^dl{rv)) + Ve^5|(ru) 

- rdl{^)dl{rv) - 2rdR{^)d%{rv) = G^ - duj{^ru%)v 


1 


u'- 


+ -7T dR{ruRR) = (3.86 


1) - (3.86 


9), 


and so we must obtain a-priori estimates for the equation: 


(3.86) 


1 


1 


-dRR{rv^) + -r.ru^jiVuj + rdji{-.dji{rv)) = fR + g. 




(3.87) 


Lemma 3.6. There exists a unique solution v to Equation (3.87) on the domain (O,0o) 
{Rq, Rq -b N) subject to the boundary conditions v,vr = 0 at {i? = Rq, Rq + N} and the 
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initial condition ?; = uq at {w = 0}. This solution v satisfies the following estimate, uniform 
in N: 

l<9R(rna;)P+ sup J r\dRR{rv)\‘^ < j J f + J j 9^{R-Rof (3.88) 

r\dRR{rv)\^. 


+ 


a ;=0 

N 


fiM i-N ^0 

Proof. Define an inner product by [[u,?/;]] := / dR{rv)dR{rw) + / —^rvw. All of 

J Ro J Ro ^ 

the properties of inner-product follow from the properties of the integral, aside from non 
degeneracy. Supposing [[t',u]] = 0, by the positivity estimate, (1.35), J \dR{rv)\‘^ = 0 
|9j?(r'u)| = 0, coupled with the fact that u = ur = 0 at {R = Rq, Rq + N} implies u = 0. 
Let Bj represent an orthonormal basis for with respect to this inner product. The weak 
formulation of (3.87) reads: 

[[vu,ej]] + j ^dRR{rv)dRR{rej) = - j fejR + j gej. (3.89) 

k 

Define R) = E aik{oj)ei{R) . Inserting this into the weak formulation above enables 


2=1 


us to solve the corresponding ODE for the coefficients aik ■ Multiplying the weak formulation 
by di^Oik and summing over i yields: 

[[vtvt]] + j ^dRR{rv^)dRR{rv^) = - j fv^R + j gv^. (3.90) 

In order to pass to the limit as /c —)• oo, we must obtain a-priori estimates for the above 
equation. We relabel by v for this purpose. In order to obtain the desired a-priori 
estimate, we multiply equation (3.87) by rv^j and integrate by parts: 


j -dRRirvc^)rVoj + ^v'^Urrv'^ > j \dR{rvuj)\" 


(3.91) 


by the Positivity estimate, (1.35). Next, we have 




J:^\dRR{rv)\ 


Ro 


-.dRR{rv)dR{rv^). (3.92) 


On the right-hand-side, we have: 

J fRrVu; + J grvuj < J + J g^{R-R(}f + ^ J \dR{rVuj)\‘^. (3.93) 
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Using Gronwall’s Inequality and integrating yields; 


poo 

+ sup / r\dRR{rv)\ 
JRo 


< 


f + I I 9 {R-Ro) + / r\dRR{rv) 

J LJ=0 


(3.94) 

□ 


Next, the following high regularity estimate is established: 

Lemma 3.7. 


+ sup J r\dRR{rVuj)\‘^ < J J + J J 9^{R-R(}f + r\dRR{rv{0, ■))] 

p p poo 

/"+ / / / r|5,jR(ru,(0,-))p. (3.95) 

J J JRn 


Proof. Differentiating the equation (3.87) once in io yields: 


-dRR{rv^uf) + ( ) rvui+ [ ] rvu,u> + tPrr ( ) dRR{rv) 


u'^ 


(3.96) 


+ rdRR -^dRRirVuj) ] = Ir^ + Qu- 


u'^ 


In order to obtain a-priori estimates of the equation (3.96), we multiply the above equation 


by remaining cognizant of the boundary condition that vr = 0 ^ Vi^^jR = v^^R = 0 at 

{R = Rq, Rq -|- A^}: 

- dRR{rVujuj) + dRR{rv)^ (3.97) 


+ rdRR ( -^dRRirv^) 
VR 


rvux. = 


Jruj + 9u} 


rv,^ 


First, applying the positivity estimate (1.35), 


- J dRR{rVujoj)rv^uj + > J \dR{rv^^)f. (3.98) 


For the second term in (3.96), we have: 


^RR 


r Vu,Vi.u, < 


^ ( f R-^vl^ ' < 


|9h(?^<.,;)|' 


1/2 


\dR{Vu,oj)\‘' 
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(3.99) 


<N{6) / + 5 / \dR{vojuj)\‘^. 


For the fourth term in (3.96), we have: 


rdRR dRR{rv)^ = 2y/e j 


n ; dRR{rv)dR{rVu. 

lU t 


+ y dRR{rv)dRR{rv^^) = (3.100 


1) + (3.100 


2 ). 


(3.100) 


For the fifth term in (3.96), we have: 

j Qrr (^-^dRR{rVuj)^ r'^Vajco = j ■^dRR{rv^)dRR{r‘^v^uj) 

= j -^dRRirVuj) {2yfedR{rv^uj) + rdRR{rVa,oj)) 

= y ^ \dRR{rVuj)\^ - j \dRR{rvu;)\‘^ + y ^dRR{rvuj){2y/edR{rvujco) (3.101) 

= dslolj l) + ( [3lMj 2) + ( |3lM| 3). 

Finally, for the right-hand side, 

y fihjrvu;u; + y gturvujuj < j + j gl{R - Rof + S y dR{rvuiujf. (3.102) 

( 3.101[ 3) and (3.100 1) are estimated through Young’s inequality. For (3.100 2) , we write: 

dRR{rv)dRR{rVujuj) = [dRR{rv)dRR{rVuj)] - dRR{rVuj). 

Applying Gronwall then yields: 

y y |9i?(ru^c^)p+sup y r\dRR{rv^)\^ < j j f+j j g‘^{R-Rof+j r-|aRR(ru(0, •))P+ 

/ [ f^+ [ [ gl{R-Rof+ [ r\dRR{rv^{0, ■))\'^ + f dRR{rv)dRR{rv^) 

J J J J J TIq j u)=9q 

- / dRR{rv)dRR{rv^). 


I Ld =0 


The hnal two terms are estimated through Young’s inequality, where we recall Lemma 3.6 
\dRR{rv)\‘^. 


to estimate N 


' llj= 0 , 9 o 


□ 


We now obtain the weighted variants of the above two lemmas. Notationally, depict the 
weight by Pm{R) = {R- ^o)™- 
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Lemma 3.8. 


Pm\dR{rVuj)\‘^ + sup jPmr\dRR{rv)\^ < j jPmf 


+ 


g^Pm+ 3 + / rpm\dRR{rv)\" 

J aj=0 

(3.103) 


and 


+ 


Pm\du,dR{rVui)\^ + sup jPmr\dRR{rv^)\‘^ < j jpmf 

/ / glPm+3+ / rpm\dRR{rVuj)\‘^ + / rpm\dRR{rv)\ 

J J J a;=0 J u)=0 


g Pm+3 + 


Pmfu) 

(3.104) 


We will proceed in several steps to establish Lemma 3.8 As these are a-priori estimates 
and we eventually plan to send N —)• oo, we work in the domain R G {Rq,oo). Dehne 
w{R) = {R — Rq)'^ on + M] and w{R) = M™' for R > M + Rq. Note also that 

w{Rq) = 0, which eliminates boundary contributions from {r = i?o}- 

Claim 3.9. sup J w{R)\dRR{rv)\‘^+J J rw{R)\dji{rv)\‘^ < J Pm\dRR{rv)\‘^+ J J fpm+ 


g Pm 


Proof. Multiplying equation (3.87) by —w{R)dRR{rv) yields: 

wdRR{rv) = j {Jr + g)wdRR{rv). (3.105) 


T'li 1 

dRR{rv^) -- rd\{^d\{rv)) 


Integrating (3.105) by parts yields: 

dR{wr 


dc J wlORRirv)]"^ + J ^^^^^dUrv)d]i{rv) + J ^\dUrv)f < J {f+g^)w + J, 

(3.106) 


where J contains: 


j = - r-^ 


VujwdRR{rv) + Jr ^|5|(ru)||5|(ru)|-h Jr "|9i?i?(ru)p. (3.107) 


Since dR{rw) = {m + l)^/ew -|- Rom{R — Rq) 


m—1. 


dRiwr] 


(m + l)^u; 3 


{R - 72o) 


m—1 




dR{rv)d\{rv) = / 1- d\{rv)d\{rv) + itiRq / -5^(ru)a^(ru) 
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\/e J ^l'9|(™)l^ + \/e J ^\9R{rv)\‘^ + j 




ISlS/" m2 I f -^o)™' ^ip, 2 /- m2 

+ / --n-|5 r(™)I ■ 




(3.108) 


The first two terms above can be absorbed into (3.106), and the third and fourth terms 
above have been estimated inductively. Applying Gronwall and integrating yields the 
desired lemma. 

□ 


Claim 3.10. sup j w\dRR{rvc^)\‘^ + j j rw\dR{rvc^)\‘^ < j Pm\dRR{rv)\‘^ + j j fpm 

9 Pm T / Pm\dRR{TVi^')\ I I ftjPm T 


+ 


/ U!=0 


9ujPm- 


Proof. We apply the multiplier —w{R)dRR{rVi^) to the differentiated equation (3.96): 

/ (Equation 


3.96) X {-w{R)dRR{rv^)). 


(3.109) 


On the left-hand-side, we have: 


rw , 


dR{rw) 


duj J w\dRR{rvu;)\''‘ + J ^|9^(™a;)|^ + J dR{rvui)dR{rVaj) + J + 1. (3.110) 

Here J = y dRR{rv)^ dR{rwdRR{rvuj))+j dRR{rVui)dR{rwdRR{rvuj)) 

and d < / v‘^ujR + l<9RR(™aj)P- J and the third term of (3.110) can be treated through 
Young’s inequality and induction after noticing that dR{rw) = {m+l)^/ew+Rom{R—Ro)^~^ 
as in the previous lemma. The desired result now follows from Gronwall. 

□ 


Claim 3.11. 


9 Pm+3- 


t(^|9H(™a;)P+ sup / wr\dR{rv)\‘^ < / rpm\dRR{rv)\‘^ + / / fpm 
J J = J J 


+ 


Proof. By Claim 3.10 the quantity Or {^dRR{rv)) is integrable and so integrating equation 
(3.87) from i? to oo to yield: 

POO rpn.O / T* \ 1 poo 

-dRirVuj) + ~ ^ Jr 

Multiplying the left-hand side of 3.111 by —w{R)dR{rVuf) gives: 


POO ^/lyO / T \ 1 

dR{rv^) - / - dR [^d\{rv) \ + \/i^ 5 |(rn) 

J LL X ci / LL 


wdR{rVuj) 
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= J w\dR{rVuj)\^ + J ^d‘ji{rv)w{R)d]i{rv^) + J, ( 3 . 112 ) 

where 

^ = y ■^dR{rv)w'{R)dR{rVuj) + j ^dRR{rv)wdR{rVuj) - j wdR{rvu,) J 

(3.113) 


J can be controlled by Hardy’s inequality, Young’s inequality, and induction for the 
rw' = y/emw + Rom{R — i?o)™'~^ factor in the f ■^d‘^{rv)w'{R)dR{rVu)) term, as in the 
previous lemma. The result then follows from Gronwall. 

□ 


Claim 


im3.12. y j w{R)\dR{rv^uj)\‘^+sup j w{R)dRR{rVoj) < j Pm\dRRirv)\‘^+j j fpm+ 

9 Pm T / Pm\dRR{rVijj)\ j I fioPm I I 9u)Pm+3 
J ijj=Q J J J J 


Proof. We differentiate equation (3.111) in lo to obtain: 


- dR{rVujuj) + j ~ (3-114) 


— fid + 


IR 


9id- 


Multiplying (3.114) by —w{R)dR{rVu]ui) yields on the left-hand-side: 


dRirVcdid) - duj 

Jr 


oo ^ ru^ 


RR 


r 


Vdd'^ - Br^ (^^dlirv)^ + yfedud ()y* 9 |(ru)^ 


wdR{rvuiud) 


= J w\dR{rVujuj)\ +duj J ■^w\dRR{rVuj)\ + 

= ( |3TT^ 1) ( [3TT^ 2) ( |3TT^ 3) J, 

where 


vw^ (li) 

- ^dRR{rVuj)dR{rVajid) + J (3.115) 


u'^ 


J = 

+ 




rw'{R)dRR{rv)dR{rVuiid) + \/e 


dRR{rv)wdR{rv 

OJUj) 


\/e y ^dRR{rv^)dR{rv^cd) + j wdR{rVujid) j d, 


ru 


RR, 




1 


~ \ ~n ] rw{R)dRR{rv) . 


(3.116) 
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Lastly, we use that rw'{R) = Rom{R — Ro)^~^ + ^mw{R) to estimate (3.115,3) induc¬ 
tively, as in the previous lemma. The result of the claim now follows from an application of 
Gronwall. 

□ 

Proof of Lemma. The estimates in the claims above are uniform in M, enabling us to send 
M —)■ oo to obtain: 

(i?-/Zo)”"|5R(r'(;a;)P+ sup j {R - Ro)"^r\dRR{rv)\‘^ < j j Pmf + j j g‘^Pm+3, 
rpm\dRR{rv)\‘^, (3.117) 


+ 


/ a ;=0 


and 


< 


{R - Ro)"'\d^dR{rv^)\‘^ + sup J (R - Ro)"'r\dRR{rv^)f 

+ Pmf^+ / / glPm+3+ / rpm\dRR{rVui)\‘^ + / rpm\dRR{rv)\ 

J J J J J UJ=0 J aj=0 

(3.118) 


Pmf + I I g Pm+3 
2 


Writing Pm ^ 1 + (72 — Ro)"^, h follows by pairing (3.117) - (3.118) with the unweighted 
estimates in (3.88) - (3.95), we can upgrade the weights on the left hand sides of (3.117) 


and (3.118) to pm = {R — Rq)^, thereby establishing our lemma. 


□ 


3.3.2 Boundary Estimates 

In the following theorem, we use the stream function formulation to give estimates of v and 
on the boundary {w = 0}. 

poo 

Lemma 3.13. / r"'(R — Ro)"^\d^^ {rv{0,t)) < C + ||r"''^^’''^(72 — 72o)™'^^'Mi||^fc+3, 

Jro 


poo 

and / r^{R-Ror\d’f+\rv^)\^ < C+||r"/2+i(i2-i2o)”^/'ni||^.+3+||(72-i2o)-”^nL.(0,r( 

JRo 


\m- 


Proof. Our starting point is equation (3.82). 
pR pR pR 

/ u{uj,9)d6, so that 'f’uj= Uuj = — 

J Rq j Rq j Rq 

define 


Define the stream function 'ip{oj,R) = 
dg{r{6)v)d6 = —rv, and ^pR = u. Now 


w = + u^ipR, 


(3.119) 
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so that: 


Wuj = -u^R'ip - + u^ipRuj = -u^jR^’ “ UR{-rv) + U^U + U^U^j (3.120) 

= + F - El- Eo + ruRR - r(t>° + vI)ur. 


According to the definitions of Eo,Ei,F in (3.79) and (3.81), / r"'R^d^E{0,ty + 

_ _ JRo 

EQ{0,ty + Ei{0,tydt < C, where the constant C depends on the previously constructed 


profiles. Now we estimate the boundary data of w in terms of ui using (3.119): 


2 

Rk^ 


poo poo poo poo poo 

/ w{0,tydt< / + / (u^yufdt < / Uidt-, and / < | llZi 

J Hq J Rq J Rq J Rq J Rq 

(3.121) 


where we have used the rapid decay of for Hardy’s inequality. We use (3.120) to do so 
similarly for Wi^: 


L 


Ro 

POO 

iRo 


Wu}{0,tydt <C+ \\rui\\j^2‘, and 


r^{R - Rord^w^iO,tydt <C + ||r”/^+'(i? - i?o)”^«i|| ^2. 


«/2+Hr_ 


(3.122) 

(3.123) 


Now we use (3.119) to express: y = 


w{u}, t) 


of both sides, and using Hardy yields: 


' Ro 


{u^y 


dt. Differentiating, taking the norm 


li^''(Ro,oo) 


+ \\r^/\R-Ri,r/^wU0,-) 


\Hk(Ro,oo)- 


J r^{R - Ror\d^+^d^yy < \\r^/^{R - R^r/^wy), •) 

This yields: 

poo poo 

/ r^{R- Riir\d^y^^ {rvy),t))ydt = / r^{R- Rord^jy\y^{yt)y dt (3.124) 

J Ro J Ro 

< \\r^/yR- Ror/Myy\\lk EIK^HR-Ror^^wUOr)\\l, <C + \\r^/^+\R- Ror/^ui\\l,+2. 


We now estimate the of v^j norm on the boundary {w = 0}. In particular, we start 
with: 


rv^ = ytoLo = u 


pR 


w 


' Ro 


,0'i2 


dRirVa,) = dujuj Ur 


pR 


w w 

{u^y 


r'^iR- Ry^\dR{rv^)y < Jw^E jwl + j r^iR-RnTw, 


m „,,2 

iJUJ 
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poo 

C{l + \\rui\\l,)+ r^{R-Rorwl^. 
JRo 


(3.125) 


Differentiating (3.120): 


Wu,oj = -u'i^R'ip - + F^- {Eo + Ei)^ + ruRR^ - rd^ {vp + vl) ur - r(i;° + vI)ur^ 

(3.126) 


(3.126 


= ( [3l26l l) + ( [3l2^ 2) + ( [3l2^ 3) + ( |3l2^ 4) + ( [3l2^ 5) + ( [31^ 6) + ( [3l2^ 7). 
We estimate each of these terms: 

.4) < J+ JiR-RoY^^l + C{u%v% 


1) + (3.126 


2) + (3.126 


^) = j Fl< j+ j + j + j 

+ J {ulv^f + J + J (u°t;°)2 + J + J {u^p^vlf + J {ulvl^f + J {u%J^ 

+ J {UpfiuLuf + J {PpJ^-^ 

7)< Jr^+^{R-RorR%R. + Jr^+^iR-Rorul + r^+^iR-Rorulj,. 


5) - (3.126 


All of the terms in (3.126 3) above are bounded by a constant C, using the estimate on 


vl, aside from the term. For (3.126 5) - (3.126 7), we use the divergence free condition. 
This establishes the desired result. 

□ 


3.3.3 Construction of Prandtl Layer Solutions 

We first define v{uj, R) = v{uj, R) — ~ Ro)'uluj{w, Rq). Here x is a cutoff function 

near 0. Then v{oj,Rq) = v{w,Rq) = 0, and vr{uj,Rq) = vr{uj,Rq) — -^ul^{u!,Ro) = 0. 
Since the equation (3.87) above is linear, we easily find: 


1 


- dRR{rv^) + -Rru^RRVu + tBrr -^Brr (rv) = fR + g, 




1 


(3.127) 


where / and g are defined: 

/ = duj{\)ruRR - ^/e\^]i{rv) + 2rd]i{rv)dR{\) + ^9^ (r(u° + vI)ur) (3.128) 

'U U 'LL 'LL 
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+ 



URr{vl + vl) 


(/•I) - (/-S), 


9 = -du;{Xru%R)v - dR^{^)ruRR - ^/e^R{^)^Urv) - 

LL LL LL Uj 

- Qr da; {r{vl + vI)ur) - dojR {uRr{v^ + vl)) 

+ daj{^){FR — Eir — Eqr) + ^ {Eroj — EiRoj — Eqr^) — da; {^) {u^rU + U^Ur) 

- ^ {da; {u^rU + U^Ur))) - gfifiCr x(-^ - Ro)ula;aj{^^, Ro)) 

+ XUeujco + i^dRR {^^dRR I r ^ xUeojj j ■ (3.129) 


We define |||f llP = 



(3.88) 


|i;||P< 


+ sup [ —prr|9/jij(ru)P. According to estimate 

[o,eo]J 

+ (3.130) 


Using the definition of v, we record the following consequence of the divergence free 
condition: 


2 ^ -2 I f R Ro \ 2/ 1 \2 I 2-2 2 2/ 1 32 


< < ev^ + 


Ro 


X^(nL)"+r^4 + ^Vfer + 


R Ro \ /\ 2 / 1 \2 


Ro 


birWe^r. (3.131) 


+ dR (wr)^ + dR (^rxuia + R^^ ~ ^o)'«L(^i 


(3.132) 


We must now give estimates on the terms in / in terms of |||u||| in order to apply the 
contraction mapping principle. First, we relate u and v to |||u|||: 


'^RR ^ 


< 


'^RR + 


J dRR < 0 osup J bP + ||?xLlli 2 (o, 0 o) 


< 


^111^ + II^LIIl2, 


eu^ + e||Mg^^||^2 < 6*0 / / + ^II^Llli^ — ^olll^IlP + ^ll'“Llli2, 


(3.133) 


(3.134) 


Ro) 
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ui < 


/ / ^i_R + ^0 / lufjji<C+ I I \dRR{v)\'^ + \\u^^\\'j ^2 < C + 9o\\\v\\\'^ + \\u^^\\'j^ 2 , 

(3.135) 

2 < ll„,l l|2 , / / ^ \\„.l ||2 , dm2 


u: 


r r /*oo 

«Llli2+/ / |5i?(rl;)|2 < +00 / |9^(rv)(0,i?)| 

J J JRo 

+ ^o// ^ II^^Lllia +^o|||w|||^ 

< J Ui + 00 J J '^‘co which has been estimated above; 

KrI"^ < [ [\dR{rVa,)\‘^ + \\ul\\l2, 

\dR{rVaj)\‘^ < \\ulJ\l 2 + [ [\dR{rVaj)\‘^. 


u, 


(3.136) 

(3.137) 

(3.138) 

(3.139) 


Now we turn to the terms in /, which are given in (3.128): 


(f.l) For ruRR we use equation (3.82): 

J J \dui^)\r‘^\uRR\‘^ 


< 


\u 


/O 

|2m,0 


012 I 


+ u^Uoj + ^/eu^v + r(Up + vI)ur + ru%v + Ei + Eq - 


:= (/.l.l)-(/.1.8). 

Each term in this equation is bounded by C + C'6lo|||u|||: 


(f.1.1) : / I R-^u^< I I R-^[u{0,Ry + eo - ^ l^il + + 



(f.l.2) : J Jr < J J^ + + + 

II^LllL2(o,eo) + ^0 J \''^Rr\‘^^ 

(f.1.3) : 




R-^ev^ < 


R-'^ev^+ / / i^-"e|uLl^ 


(f.1.4) : J J{v^, + vlful< J jul<j J uIr + Oo 

<C + eo(^J J\dRR{rv)\^ + J 


u: 


ujR 


(f.1.5) : J J R-^v^ < 


R-^v^ + / \ul\\ 
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(f.1.6) - (f.1.8) : The forcing terms decay rapidly and therefore < 


R-^\F -Ei-Eo\^ < C. 


(f-2) e j j\^\‘^\dRirv)\‘^ < e\\\v\\\‘^ + e j j {Orr XulJ^ |^ 

(f-3) j j r^\dRirv)\‘^\dRi^)\‘^ < j j R~'"\dRRirv)\‘^+j j R-'^IBrr P 


(f-4) J J\^\'^\da;{r{v^ + vl)uR)\‘^ <\\vi 


ecj M oo 


Ur, 


(f.5) / l\du;{^)\‘^ujir^{v° + vlf< 


R-^u%\vy + I I R-^ul{vlf 


Ur- 


Combining (/.I) - (/.5) with (|3.133|) shows that j j f < C + 6»o|||i;|||^ + II^^Llli 2 (^=KQ) 


We now give estimates on f f g‘^{R — Ro)^, recalling the definition (3.129): 


(g-1) J {u"^ + u'rr + \dRR{rv)\^ + uIj + uji ++ uj^) < C+9o\\\v\\\‘^+\\ul^\\‘j^2(^j.=Ro)’ 

(^) +^e)l^ ^ J J ^R- 

Here, we have used: f f R~'^lvg^l^uj^ < \\R~‘^vl^\\oo f f uj^, and ||i?“”Ue^||oo < 

||i?-X'lb3 <e-||r-ui||^3 <C. 

(g.3) I l\dU^)\^\FR-E,R-EoR\^ +11 \^\^\Er^ - E,r^ - Eor^\^ < C hy Enlei 

bounds, 

x{R - RoWe^^iuj, Ro)^)\‘^: 

j ul^^iuJ,Rof < j j driul^j"^ = j = j j {vl^ + rvl^^) {vl^r + drirvl^J) 

< ||r"'Ug||w2.p||p”?^ellvF3>'j ^ for K arbitrarily small. 


The boundary terms in (3.129) can be estimated by 


r-^o 


.1 |2 




+ 


17/* 1^ < 

I Ueujuj I ^ £ 
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Finally, it remains to give estimates on 


' Ro 


r—ndji{rv{0,R)) from (3.130) 


pCX) 

/ r 

Iro 


1 poo poo 

-Qd]i{rv{0,R)f < / r\dRR{rv)\‘^ + / r\dRR 
^ jRn JRn 


R -^0 \ |2| 1 12 

X ) I WeJ 


iRo JRo V -^0 

< C + I1^4 + < C* + ||Fi||j^4 


Therefore, we can close the contraction mapping argument. Because the estimates are 
uniform in N, we can let —)> oo to obtain a global in R solution. We can repeat this same 
argument with the weights of R^. Indeed, the only terms above sensitive to a weight are 
contained in /, and we treat them here: 


R^{v^p + vlfu%< j I R^u%<C+ I j R^\dRR 

^\^R!^\dRR{rv)\^ <e j j R^\dRR{rv)\‘^ + e j 


rv\‘^ + / kLl^ 


-^\‘^R'^\d^{r{v^ + vl)uR)\ ^ I I li- aji 


< 

r\_/ 


R^ul + 


pm 2 
^ ^Rui- 


(3.140) 

(3.141) 

(3.142) 


Again through contraction mapping, this establishes: 

Lemma 3.14. J J ~ Ro)'^\dR{rVa})\‘^ + sup J ~ ^o)^'i"\dRRirv)\‘^ < e~'^ for k > 0, 
arbitrarily small. 


Uniform estimates are obtained via the calculation: 


2 i-e 

;r. \2 , / IQ 


u{uj,Rf = \ui{R)\‘^ + dR{rv)j < \ui\'^ + J \dR{rv)\ 

/•do rR 

<C+ / dRR{rv)dR{rv) 

Jo J i?o 

<C + [ f R-^\dR{rv)\^ + [ f R^\dRR{rv)\^. 


(3.143) 


Similarly, using the fact that v{uj,Ro) = 0, we have 

rR 


pit roo poo 

v{uj,R)‘^= / vvr < / R~"^\v\‘^dR+ / R^\vr\ 

J Rq j Ro j Ro 


<1 R-^R^\\vrr\\12+C+ I I \vrJ\ 

' Ro Jo J Ro 


Ro 

6o poo 


(3.144) 
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This yields: ||f||oo + Hi^lloo < Ce We now obtain higher-regularity estimates. The 


starting point is (3.95), and as such we define 

j j \dR{rv UJUJ )p + sup j r\dRR{rVaj)\^ 


hlll' = 


(3.145) 


Since / and g have already been estimated, we compute and Quj: 

fuj = dujuj{X)ruRR + du;{\)ruRR^ - Ve^dj^irvuj) + 2rd]i{rVuj)dR{\) 

+ 2rd\{rv)duiR{^) + {r{Vp + vI)ur) + URr{vl + vl) 

+ u^Rr{vl + vl) + URvd^iVp + vl) + d^{r{Vp + vI)ur). 

(3.146) 


9uj = -d^{-f^ru%ji)v - duji-^ru%ji)vuj - dR^{-^)ru^RR - dRaji^{-^)ruRR 

- VedujR{^)dUrv) - VedR{^)dl{rv^) - 2rd]i{rv^)dl{^) 

- 2rdji{rv)dRRuj{\) - d^R ( ^ (r(u° + vI)ur) - Sr f ^ ] dojcj {r{vp + vI)ur) 

U \ 'LL J \ U J 

- dujujR {uRr{v^ + vl)) - d^R {uRr{v^ + vl)) 

+ duj (d^{—^){FR — Eir — Eqr) H— {Fr^ — Eirj^ — Eqr^) 

\ VR vr 

- duiuj{^) + uI,ur) - 2d^{^)d^ [uIjrU + v^ur) - ^ {uZrU + u^ur))) 

- dRR{r ^ - Ro)ul^c,Loi<^^ Ro)) + duj ^ Xui,uj 

+ + rd^RR (J^p^RR ) • (3.147) 


We now treat the terms in (3.146) 


f < 

J UJ — 


+ e 


R {v‘rR + VRRuj + \9RR{rVu]\‘^ + uji + ujii^) + J J \^\‘^\du}u}ivluR)\" 


u'^ 


(3.148) 


All of these terms can be estimated in terms of the norm, with urr^ being estimated by 
taking of Equation (3.82). We now address bearing in mind that all of these terms 
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are accompanied by rapid decay: 


+ 


9t{R — Ro)^^ J J ^ ^ 


R-" \vLful + 


dui ( dui{—^){FR — Eir — Eqr) H—^ {Eruj — Eir^ — Eqrui) 

' vy 


+ u^ 
(3.149) 


All of the above can be estimated in terms of |||w||r. The most delicate boundary terms 
< C{dQ)e~‘^, and the most delicate interior term is: 


J r=Ro 


-n„,l ||2 

eujuj \ \ oo 


R~'"u% < \ \R-^vl\\l 


eWm 


R--u% 


<e-\\r^vl\\%. 


Ur- 


(3.150) 


The latter boundary term in estimate (3.104) has been estimated before. Therefore, we 
must estimate the Vt^ boundary term: 


[ r{R-Ro)^\dRR{rv^)\‘^ = [ r{R - Ro)^\dRR (rd^ (v + 
J u=0 J u)=0 \ \ 


ito 


< / r{R-R,r\dRR{Tv^)\^+ [ rR^\dRRi^^-^xul^JO,Ro))\^ 

Juj=0 Juj=0 -f^O 

< C + ||r3/2iiW2^^|||3 + ||i?-™uLjO,r(.))||Hi + l«LjO,i?o)P 


< C(0o)e-'/^ 


(3.151) 


where for the final inequality we use the same calculation as in [GN14t page 25]. Therefore, 
we can close the contraction mapping argument, and again let —>■ oo yielding: 

Lemma 3.15. j j{R - i?o)™|9i?(rUa;^)|^ + sup j(R - Ro)'^r\dRR{rVuj)\‘^ < C'(6'o)e"^ 
where the constant C{9o) could depend poorly on small 9q. 


3.3.4 Cutoff Prandtl Layers 

We define our Prandtl-1 layers (uj,, Vp) by cutting off the previously constructed layers 
{u,v): 


rR 

Up{u,R) :=x{V^{R-Ro))u + y/exW^{R-Ro)) u, (3.152) 

J Ro 
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R) := x(\/e(^ - Ro))v. 


Here, x is a standard cutoff function which equals 1 on [0,1]. (Up, Vp) satisfy the divergence 
free condition: du,Up + dji(rVp) = 0. We also have the following estimate for Upi 

rR 

,,1| ^ . 77/' r? 

j-p 11 oo ^ 


\^ex\V~e{R-Ro)) f nl\<V~e{R-Ro)x'H" " 

J Rct 


(3.153) 


We have the following lower-order estimates of ul and v, 


1 . 


r^{vUY < e 


- 1/2 


R^ivL^RY < 6 


-5/2 


r^ivLY < e-'/" / / R^{vLrY < 


rYuU" 


< e-i/2 / / R^{vUrY < 


pRR) 


(3.154) 

(3.155) 

(3.156) 


\dRirVp)\ = 



dRirvYdRuirvl) 


R 


< £-1/2 


< £-1/2 


J j \dR{rv)\\dRR{rv)\ 

11 R-^\dR{rv)\^ + 11 R^dRRirvY 


2 \ < ,-1/2-h 


(3.157) 




prY ^ j j RpujRi^: R)\‘^ < C + J J \u^rY — C + j j \dRR{rVp)\‘^ 

< €-^, (3.158) 

(3.159) 

(3.160) 






Wlu? < e 




Ro~^~7e 

I 1|2 


Vpf ^ \\Vp\\l- 


Ro + tTJ 


•1"2 I I dRdco <e-" 2 -^ 

^ Rq j j Rq 

The weight of r” can be added in because on the support of Vp and Up, r is bounded 

uniformly. For the same reason, these weights can be added in for the uniform estimates. 

We summarize the results of the Prandtl-1 layer construction in the following: 

Theorem 3.16. Let Up,Vp denote the cutoff Prandtl layers described above. Then 


R^\dR{f’v^)\'^ + sup J R^r\dRR{rVp)\‘^ < C{6o)e ^ for k > 0 arbitrarily small, 

J J ii™|9pj(ri;i^^)P-hsup J R”^\dRR{rvl,J\‘^ < C'(6»o)e-2, and | ui| |oo < C'(6'o)e-'^, 

where C{6o) is a constant which could depend poorly on small Oq. 
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3.4 Estimates on Profile Errors 
3.4.1 Angular Errors, R^\ 


There are three lower order error terms from the profile constructions, which come from 


(3.22) for the Prandtl-0 layer and (3.72) for the Euler-1 layer and for the Prandtl-1 layer by 


inserting (3.152) into the equation (3.78): 


poo 

\e 2 ) + J Eb{uj,6)d9) + (Prandtl-1 contribution) = (3.161 


1) -f (3.161 


2) -b (3.161 


3 ). 


(3.161) 


First, we have 


(3.161 1) = ||r™e 2 ||L 2 < \\r^u%R'^\\L 2 \\d^irv^)\\L 2 + e\\r^ul^R'^\\L 2 \\u]l^^\\oc < Ce 


(3.162) 


Next, we have: 


Eb{io,e)d9ydRdu = e -^2 J J Eb{uJ,9)d9y drduj 


<6-5 / x(^)2da;<e5. 


(3.163) 


To estimate (3.161 3), we write: 


Prandtl-1 angular error contribution =: 

f-R \ / rR 

0 / 


Xul -b y/ex' / ui -bn° xuluj + Vey' 


+f{vl+vl) + 2v/ex'4 + ex" J ^ 

- r {xuIrr + ex"u^ + 2\/ex'u^pj -b j -b 2ex" + V^XuIr^ -{E-Ei- 
Eo) 

= {l-x)iF-Eo-Ei)+y/^x' (ul J u^ + u^ J ul^ + 2(u° -b u^)u^ + -bex" (^(r^p + vl) J 

/ rR \ 


u, 


■p I • 


We give estimates on each of the terms in the above expression: 

r^{l-xf\F-Eo-Ei\^<d^, (3.164) 

nx'{V^-)?H + vl)ul + ulp, + ul I |x'(Vi^)pi?-" < 6i-^ 

(3.165) 


+ 


R 

Up -b 3Up 
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el I \x'{Vi-)Wf\ J y (3.166) 


ix" ((f? + vl)J + e^/V" (y uj' 

<«"//(x"(^-))"(y^«;,) •'’"<e2-«yy(x''(^-))"flV 


X 




(3.167) 


Now, we must estimate the higher order terms from the expansions (3.1) - (3.7) : 


e^-order error: 


^™ke+^plVL + ^pa;)P < ll^e+^pIlL 


nv^, + vl\^\ul\^<\K{v^, + vl 




\u^ P < e-1/2 
I ^er I rv.^ ^ 5 




l^er + '‘^prI ^ ^ 


1 

v/i ^ 


-p/ rpi j j (“? + “)l)'^Se “ 

We dehne: 

52 _ / /'„.o,.o , „,0 d ^„,o I „,l^^ I /„,0 I „,l^„,0 


(3.168) 

(3.169) 

(3.170) 

(3.171) 

(3.172) 

(3.173) 


poc -p 

Jr ^ ^ + vt)) + {vp + Ug^)u^R - - vluRdt 


r’Ko+y? o /■^o+“^ 1 

—^uiupdt — / ——uiu^dt. 

r{t) P P Jr r{t) P ^ 


JR 'K‘') JR 

Using the rapid decay of Prandtl-0 prohles, after taking d^j, the first integral is bounded 
by 7?“"' for arbitrarily large n. We must therefore treat the d^j of second and third integrals: 


cRo + tt; 


Ai+ T""*—+ U"* 1 ,.1 

r(t) P‘^ P ^ r{t) P P^ ^ 


0 




JR 'lU JR JR 

The middle term above is bounded by For the first and third terms, we use the 

divergence free condition and integrate by parts, bearing in mind that = 0 due to 
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the cutoff: 


iR 


^0+3; 2 f 

—= I - / 

JR 


Ro+^ 


r{t) 


R 

< e-'Kl + 


r(t)vl8,(^ul)dt + tijujl 


POO 

JR 




Ji? ' Ji? 

We place the terms above into integrals, the most delicate being; 




r{t)vpdt{—^)dt j dRduj 


i?o / 


-Ro+^ / rKo+3r 


Ve| 


i?o -Ro 

Rq~\~~^ f /“^o+l 


r{t) r(t)2 

2 




-Ro 


a -Ko + 1 7/^ \ / 

\uerir)\ + \^\drj dRduj < 


^^L ,0 |2 I |„.0 |2 I I 1 |2 1 I 0 |2 I I 0|2 1/2 


^r2+'5di?<||u0"2 


— M *^6 M OO 


(3.174) 

hs‘^R<^-'''^ r (3.175) 


We emphasize that estimate (3.175) is the most delicate and requires the weight parameter 


5 of to be strictly less than 1. 

e^/^-order error: 

r-|ui|VL + + I I r^Wl +< + < + vlJ" < 


T^\u^ -L 7/^ -I- lA A 

\ ^err ' ^er ' ^ecjcu \ 


puj ' ^euj\ 


\m 


\m 


r^\u^ |2 

' I ^pUJiO I 


11 \dRirv^)\‘^ < 6-3 


e^-order error: 


mi l |2 < -1/2-h 

' I ^puj \ ~ 


-«:-l/2 

(3.176) 

< 6-3/2, 

(3.177) 


(3.178) 


(3.179) 
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3.4.2 Radial Errors, R^\ 

The higher-order radial contributions must be estimated: 


gi/2 order error: 








+ [ [ {ulv^co + + vl))^ + I f r^\vp\‘^\vpji\‘^+ < € ^ 


(3.180) 


r- (\vL\^\4 + 4\^ + K\W + K?Kr\" + + \vl?KRf 


^puj I 


< g-«-l/2 


m /I 0 |2 , I 0 |2 , I 1 |2 

r iiu„dI + \Up^\ +\Vpjiji 


)< 


^pR 

e order error 

r'^lvlrr + vlr + < (^~^^‘^\\vI\\h2 , 




„m| 1 |2 

^pULO \ — 


We have established the main result of this section, Theorem 1.10 


(3.181) 

(3.182) 

(3.183) 


r""\ul>v4 + ulv4 + v^vl^ + + vIr + + Vpl"^ < e (3.184) 

g3/2 order error: 


(3.185) 


4 Energy Estimates 


In this section, the energy estimates in Theorem 1.11| are proven. We will need to work 
with the restricted domain il]\f = (0, 9q) x (Rq, Rq + N), and obtain estimates which are 
uniform in N. The boundary conditions u = u = 0 on {(w, R) : R = Rq + N} are enforced. 


Step I: Laplacian and Lower Order Terms 


By the divergence free condition: 



e 

T ^ UuiLil 


e 


H —dco 
r 
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(4.1) 


^3/2 

9 9 '^LO ^LoR^ 

^ fj^ 

n , , 

0 = -Vrr + Vrr = -Vrr + Or - V - 

\ r r 

y/^ 6 RlijR I /—Ruj 

= - vr + -h Ve^. 


(4.2) 


In this step the terms in (1.25 - 1.26) which involve {u,v) but do not depend on the 


profiles, {us,Vs), are treated. After adding in (4.1) - (4.2) these terms are summarized; 

(4.3) 


■\/ 6 26 6 ^ S/2 ^ 

— '^RR - Ur -9tta;w + -^U -^6 ' - V^R, 

O ^ , O ^ ““-R fA A\ 

-2vrr - VR - kVujoj-\ - ^Uuj + 2^v -. (4.4) 

ry» ry>^ ry^^ ry^^ /y* 


As described in the introduction, we proceed to multiply (4.3) by and (4.4) by 


g^i+(5.y. 


2e 


y e o 0/2 e 

- URR - UR - T^UujuJ + -^U - ' Va; - VuiR 

ry* fy*^ ry*^ ry^^ ry* 

O 2Ve e 3^6 , „ e u^r 

- 2vrr - vr - 4 - + 2^v - 

ry* ry^^ ry*^ ry^^ ly* 


+ 


X 


X er^+'^u. 


(4.5) 


Integrating by parts the terms in (|4.5|) yields: 

,l+<5 _ 


I ye 2e eu 

{urr - y^R -+ ^)ur ^ = 

ry* ry» ^ ry^ ^ 


- 2e 


/ r ^~^^uuuj + Ji, 

J u)=6q 


n|r^+'^ + 2u^r“^+'^ + eur“^+'^ (4.6) 


i-‘2vRR - 


2^/evR ev 

I 


f^ + 2^)6uri+'5 = 

j^Z rpZ 


2eu|r^+'^ + e\^r“^+‘^ (4.7) 


+ 2e^uV"^+^ - [ e'^r-^+^vv^ + J 2 , 

J u)=6q 


-e I I u^Rvr^ -e / VujRur^ = 2e 


URVujr 


5+ 5 , 3/2 




(4.8) 


/ euRvr = J3 - / euRvr^ 
/ uj=Oo J (jJ= 9 q 


59 






















where Ji = 5\/e 
Putting the remaining terms into J 4 yields: 


uuRr^, J 2 = / / vvRr\ and J 3 = 2e / / URVi^r^+5e^^'^ ^ 


J 4 = —3e2 


VujUr ^■'■'^ + 36 2 


VUojT 


-1+5 


Letting J = Ji + J 2 + J 3 + J 4 , it is easy to see that: 

\J\<C{e,9o)\\v\\l, C(0o,e)~O(0o,\/i). 


(4.9) 


(4.10) 


Using the stress-free boundary condition, (1.32), the boundary term from (4.7) cancels 


with that of (4.8): 


g2^-l+5^ - 


euRvr^ = 0 . 


luj=B(} 


luj=0O 


The only remaining boundary contribution is then that of (4.6): 

di = - 2 e f 

J uj=0q 


(4.11) 


(4.12) 


The remaining interior terms from (4.6 - 4.7) are then: 

I^ = I f -h 2ulr~^^^ + €ur~^^^ + 2ei;|r^+^ + + 2e^v‘^r~^^^. (4.13) 


Summarizing the calculations from this step, we have: 


r^+^u X (4.3) + er^+^u x (4.4) = U + /3i + J, 


(4.14) 


where di is defined in (4.12), J is estimated in (4.10), and the interior terms Ii are defined 


in (4.13) 


Step II: Profile Terms 


In this step, we treat the terms from (1.25 - 1.27) which contain Us,Vs- For clarity, we 
display these terms here: 


11 ^/e ^/€ 

-UsUoj + -UsojU + UsRV + VgUR H- VgU H- UgV, 

r r r r 

11 2 1 

-UgVuj + -VsujU + VgVR -h VsRV - ^UgU. 

r r r y/e 


(4.15) 

(4.16) 
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Applying our multiplier to (4.15 - 4.16) yields: 


X 


1 1 \/£ \/e 

-UgUaj + -UsujU + UsRV + VgUR H- VgU H- UgV 

fy* y* 


UgUijjUr^ + UgujU^r^ + UsRUvr^~^^ + VgURur^^^ + y/evgr^u^ + yfeugr^uv 


X 


11 2 1 

-UgV^ + -Vg^u + VgVR + VgRV - j^UgU 

r r r Je 


eugr^vvuj + evgcjv^uv + eVgVRvr^~^^ + eVgRv‘^r^~^^ — 2^/tUgUvr^ 

< C(0o,e)||r;||| where C(0o,e) ~ O(0o, Ve). (4.18) 

We have used the following uniform bounds on the profiles, recalling that Us{u},R) = 
ttg(r) + Up{uj,R) + ^/eul{uJ,r), and Vg{uj,R) = Vp{uj,R) + vl{u},r). 


M oo 


<c, 


||oo — W^P 
roo 


< \\v%r'^\\oc + Wvlr^Woo 


<C + \\vlr'^\\H2 <C, 


sup / Ugj^r"'{R — Rq) < C as shown below in (4.27), 
wG[O,0o] 'J 

||i^sRr"||oo < Wr'^VpjiWoc + \/i||r""yg^||oo <c + \/i||r"i;g||jy3 < C, 

11 I loo < <^(^ 0 ) by weighted Euler bounds in estimate (3.74), 
Wugcjr^-Woc < lkpa;^''lloo + \/e|loo < C* + \/e|^'''9,.(ru)11oo, 

<C + \/e||r”Ug||oo + \/e||r”Ue^||oo < C, 

ll^spr^lloo < Ve|k”'Werlloo + Hupr^Hoo + e||r”Ug^||oo <c + e||r”Ug 




< C. 


(4.19) 

(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 


Of these, Ug{uj, R) = Ug{r) + Up{io, R) + ^/eul{u], r) is the term which cannot absorb factors 
of r due to the outer Euler flow being bounded below away from zero by assumption. 


For this reason the estimate for y/e f f UgUvr^ in (4.17) is the most sensitive to the weight 
and forces the loss of one factor of r in the energy estimate. 


The estimate ||ng^||oo ^ Ikellp^ from (4.26) is given below: 

2 


( ^uj \ 2 puj pOq 

j dr{rv){e,r)de\ ^ [ul^{uj,r)) < Go J \drr{rv){e,r)\‘^de < Oq J \drr{rv)\‘^de, 
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rSo 

=► sup |Ug^(a;,r)p <9o \drrirv)\‘^{9,r)d9 := 6 »o(^(r) => ||tte^||oo < 6 * 0 sup \(pir)\. 

LO Jo r 

Since (p : [Ro,oo) —>■ M, according to the 1-dimensional Sobolev embedding: 

POO POO POO p6q poo p9o 

sup |(^| < ||vp||^i.i < / \p{r)\dr + / \drpir)\dr < / / \drr{rv)\‘^d9dr + / / \dr 

^ JRq JRo JRq Jo JRq J o 

Combining these gives 11 u].^ \ \ oo < Wr'^vlWui for a constant independent of small 9q. We 
can repeat the argument for a weight of r"". 


The profile estimates in (4.19) - (4.26) are all independent of small 9q^ with the exception 
of ||usa;||oo- Any time ||usaj||oo appears in our analysis, it must be accompanied by a power 
of e to overcome the potentially poor dependence on 9q. For the uniform bounds, this 
follows from the fact that the corresponding estimates were independent of small 9q. 

POO 

For the sup / ulj^r"'{R — Rq) estimate, we have: 


i*;e[o,0o] ■'iJo 


sup / {R — Ro)r^Ugji < C{vp, Up) + sup I {R — Ro)r'^e \u^r{r)\ dR 


[0,6o] 


' Ro 


<C + esup r^+^\ulir)\‘^dr <C+ e9Q^\\r'^+^vl\\H2. (4.27) 


I Ro 


Step III: Pressure Term 


Applying our multiplier to the two pressure terms in (1.25 - 1.27) yields: 

J+s 


Pujur + 


PRvr 


PucoT^ — {1 + 6 ) I I y/evPr^ — / j r^~^^PvR+ / Pur^ 

J uj=9q 


-<5 / / P^/evr^ + / Pur 
J J J uj=9q 


(4.28) 


Using the stress-free boundary condition at {u = which is Pr = 2euij, this boundary 

(4.29) 


contribution cancels the remaining boundary term /3i from (4.12): 


/ Pur^-2e / 

J IjJ=9q j U)=9q 

The interior term is estimated as follows: 


nn^r“^+'^ = 0 . 


\5 I / P\/~evr^\ < 6 


p2r& 


ev'^r^ ] < d9o\\v\\B\\P\\L2 


(4.30) 


{rv)\‘^drd9. 
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Step IV: Right-Hand Side 


fr^^^u<N{5) I //V+'' + 5/ lu^r^<N{5)l / /V+''+ |i;|(4.31) 


p2.„2+(5 


, 2„<5 




gr^+^v<N{5)e / / g^r^^^ + 5 / / ev^r^ < Ne / / + 0q||u 


Is- 


(4.32) 


Putting the calculations in (4.14), (4.17 - 4.18), (4.30), (4.31 - 4.32), together yields 
Theorem 11.111 


5 Positivity Estimate 


In this section, we establish the positivity estimate given in Theorem 1.13 First, we must 
prove the positivity calculation in the generality we require, that is, using the weight of r^: 


Proof of Lemma 1.12 for general S. 

r^\dR{rv)\'^ = j j r^\dR{^Us)f = j j r^\dR{^)us + Usr{^)\‘^ (5.1) 

r^\dR{^)\'^ul- j j 6r^~^{^)‘^^/eUsUsR- j jU sUsrr, 


u. 


r^\dR{rv—)\‘^ = / / r'^|nsdR(—)+ Uss( — ))P 


< 


Us 

r^\dR{-)\\l 

Us 


Us 


r^\dj,C-)\^ui+ / r^'^yuin 

J J 

(5.2) 


where we used the Fundamental Theorem of Calculus, the hypothesis that minu^ > 0, 

e terr 
\ 1/2 


and estimate (4.27) for the rapid decay of Usr. Lastly, we deal with the term; 

1/2 


r < 5-1 2 2 '^sR 
'er r V - 




2 2 2+<5 \ 

UsrV r + I 


< 


(9o||u||s ( sup f uljiiR- Ro)r'^dR] 


1/2 




1/2 


Estimates (5.1) - (5.3) immediately imply the desired positivity. 


(5.3) 


□ 
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Step I: Positive Profile Terms 




As discussed in the introduction, we apply the multiplier (r°dfi( -), —eduj( -)) to the 

_ _ Us Us 

system (1.25 - 1.27). Here we treat the three profile terms which enable us to obtain the 
necessary control of ||u||b. Explicitly, we are computing: 




UsVuj 1- 


e- I — 


(5.4) 


We now compute each of the terms in (5.4) individually: 


UsT ^^^UojdRi^) = j j UsT ^^^{-dR{rv))dR 

\/er^vdR{rv) — J J ^ J J 
^/er^v^R{rv) - J J \dR{rv)f + ^ J J r^^dR{{rv) 


For the first term in (5.5), we estimate: 


\/er^vdR{rv)\ < Oq 


,. 2„<5 


r\dR{rv)\^] <0o\\v 


\B- 


(5.5) 


(5.6) 


The second term in (5.4) is: 


vUsRr^dR{- 


Us 


6 I r^u dR{rv) 2 Wij 

vUsRT ve-hr- r v^:r 

Us Us 


VUsRT 


^-TV 2 Usr\ . ^ f f S'^sR^ /i n 2 \ ft- 


Combining (5.5) - (5.7), integrating by parts the final terms in both (5.5) and (5.7), and 


recalling estimate (5.3) yields: 


J y + r*^9 k < 6 'o|kill - y J r^\dR{rv)\'^ - (5.8) 


The third term in (5.4) is: 


—e 


r^UsVt^ duj{ — ) = -e 

u. 


r‘^u| + e / / r^^^vvuj = (|5.9 


l) + (5.9 


We estimate (^ 2): 

SUsLO 


r- -< 6»o|ksa;||oo 

Us 


evlr^ 


1/2 


evlr^ 


1/2 


< 


•^olkl 


B- 


2). (5.9) 


(5.10) 
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Putting (5.8 - 5.10) together and using the positivity estimate (1.35) yields: 


r|(9ij(ru)p + e / rvl < -(5.4) + 6q\\v 


\B- 


(5.11) 


Step II: Remaining Profile Terms 

We now apply the multiplier to the remaining profile terms and provide bounds on each 


term, keeping in mind the estimates on the profiles we proved in (4.19 - 4.26). For reference 


we include the specific profile terms from equations (1.25 - 1.27) that we treat in this step: 

(5.12) 


1 ye ye 

-UsuiU + VgUR H- VgU 4- UsV, and 

r r r 

1 2 1 

-VsujU + VsVR + VsRV - ^UsU. 

r rye 


Applying the multiplier to (5.12 - 5.13) yields: 

1 


UsujU + VsUR H- VsU 4- UgV r^dR 

r r r J \ u, 


r V 


—e 


1 2 1 ,^ 

-VsojU + VsVr + VsRV - -=UsU du. 

r rye 




Us 


We estimate each of the eight terms appearing in (5.14) - ( 5.15| ) individually: 

UsujUr^~^dR = 

<00 (I lr'\dRirv)\ 


_i+5 / ^rv dR{rv) 2 o / 1 ^ 

r ^Usuulye -hr- r vOrI —) 

Us Us Us 


1/2 


+ 


er^v^ 


1 / 2 N 


< 


(5.13) 

(5.14) 

(5.15) 


0o||r||^, where the constant C = ||nsaj||oo 4-sup / u^sr{R — Ro)r'^. (5.16) 


VsUR^dR ( - ) = 

Us 


S ( ^rv dR{rv) 2 o / 1 ^ 

VsURr ye -hr-hr vOr^ — ) 

\ Us Us Us 


<{ I I 


1/2 


1/2 

ev'^r^ ] 4- 


r^\dR{rv)\‘ 


1 / 2 ^ 


^ ^o||r||_B||u|U) where the constant is: C = Ursr’^Hoo + sup / ulRr'^{R- Ro). (5.17) 


^evsur^ ( ^ ) = 


f f Viv.r-^*‘u (VI- +rMrl+ 

J J \ Us Us Us 
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<111 eu^r^ 


1/2 


1/2 


r^ev‘^ 


1 / 2 N 


<-v/e0o||^|ls) where the constant (7 = I l^s I loo + sup / — i?o)- (5.18) 


^Usvr^ ^dRi—\ = 


(5.19 


(5.19 


(5.19 


Us 


1 ) = / / er^v^ <Oo I I er%l, 


\/e A /-fu dji(rv) o „ , 1 , 

—Wsur" ( y/e -h r - - + r^vdR{ —) 

r \ Us Us Us 


(5.19) 


= (5.19 1) + (5.19 2) + (5.19 3). 


„|5„.2 


2) = j j y/er^vdR{rv) <eo(^J j (^j J r^\dR{rv)\ 


1/2 




< do 


^) = j j Vev‘^r^^^UsdR{^) ^ (^J J 

sup j r^u^RiR - Ro)^ {^j j\dR{rv)\ 


t’^uIrTv 


1/2 


1/2 


Thus, (5.19) < We now individually compute the terms in (|5.15[): 


5 I U^jj VUsu) \ , r l\\ II ill II 1/311 I 

er Vsu}U ( ^ I ^ ye (I lusoj 11 oo T | \ usuj 11 oo) I tI 

Us ui 


(5.20) 


Here again \/e||ustj||oo is a good term despite the potentially poor dependence of Ijus 
on 6q. 


-er^^^VsVR 


-er^^^VsRV 


Vu, VUs 


2yfer 


1 

\ ^ 

ul 

( Vuj 

VUsu 

1 

1 

ul 

Vuj 

VUsuj\ 

Us 

ui J 


^ < \/e(||rUs||oo + II '^'^SUJ I I OO ) I I ^1 I s 5 

oo + 11 TUguj 11 oo ) ^011 11 _B 5 


(5.21) 

(5.22) 


< 




1/2 


er^vl 


1/2 


+ Us 




1/2 


er^2 


1/2 


< 


(5.23) 


6 'o||u||s. 

Terms (5.19) and (5.23) require precision with regards to the weight in our multipliers, 


as the prohle Us cannot absorb any factors of r. The results of this step are summarized: 

(^ + (ra <C^(0o,e)l|u||| + iV||u||i, C{6o,e)r.O{0o,V~e). (5.24) 
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Step III: Laplacian and Lower Order Terms 

In this step, we treat the terms which contain {u, v) and do not depend on the profiles 
Us, Vs in equations (1.25) - (1.26). For clarity, these terms are summarized here: 


\/e e e ‘^ 3/2 j 

URR - Ur - ^Uojuj + -^u -^ 6 '’' Voj, and 




e 2^/e e 

- VRR - VR - kVujuj H-+ -itV. 

ry* 


(5.25) 

(5.26) 


Applying the multiplier then yields: 
(Equation 


5.25|) X r^dR{^-^) — e / / (Equation 


5.27) X d^{ -) = 

Us 


^ ^ 2 3/2 \ 

-URR - -^Ur - -Z^U^UJ + -^u - ' v^]r Or[ — ) 


— e 


y/e e 2y/e e . 

-VRR - VR - ^V^uj H- ^u^ + -^v duj 


Us 
rp\-\-5 y 


(5.27) 

(5.28) 


We proceed to individually estimate all of the terms appearing in (5.27 - 5.28). 

..2, 


-urrv dfi{ -) = - 


Us 


[ [ URRV^ ( Ve— +—dR{rv)+ r^vdR{ — )\ 

J J \ Us Us Us J 


= ( [5^ 1) + ( [5^ 2) + ( [5)^ 3). 


(5.29) 


(|5.29[1) = Ve / [ urOrC^ - -) = e f f ur— + Ve f f UR—dR{rv) 

-' J J Us J J Us J J Us 

+ \/ey j URV^^^vdRi^) < ^/e\\u\\A\\v\\B, 

mv = l /“«»«(' 




= -(1 + 5)\/e UR 


r^Uuj 


Us 


URV 


1+5 UqjR 


J J r^^^URU,^dR{^) 


= ( [5)^ 2. 1 ) + ^^ 2 . 2 ) + ([5^2.3). 


^2.1) <V~e(^l I r^ul^ (^IJr^ul^ < v/i|klUI|i^l|B, 

@2.2) = «) = 4 /(i) - \ I 

(|5.29[3) = [ [ URdR{r'^^^vdR{ — )) = {2 + 6)y/e [ [ URr^^^vdR{ — ) + [ [ URVRr^^^dR{—)+ 


luj=dQ Us 
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r'^^^URvdRR . 


For the last term in (5.29 3), we use: 


sup / < C” + fc sup / r"^{R - Ro)\u\^^\^dR 

J Ro J Ro 

<C + e^ sup J r"^\ulj.j.\‘^dr 

< C + e^e^^Wr^vlWl^ <C + e0o\ (5.30) 

poo p 

sup / r^{R — Ro)\usr\'^ < C + sup / r'^{R — Ro)\ul\‘^dR < C + C{6o)e^. (5.31) 

Jro j 


term in (5.25) is: 


- 1 /* 1 

Summarizing, we have shown (|5.29|) < (5||u||^ + A^(5)||tt||^ — - / —m|.. The next 

VjO 


2 .luj=eo 



r Ur 


1/2 






er^v^ 


1/2 


1/2 


-^feuRr f^/e—+—dR{rv)+ r‘^vdR{—)] 
V Us Us J 

<%/i 

+ OoWusRr 



^l+5,,2 

r 


1/2 



\dR{rv)'\‘ 


1/2 





1/2 


< {Oo + \/e)ll^i|UII^^| 


B- 


(5.32) 


Next, 


-(^^zsUuujdR{ — ) = [ [ -^u^d^dR{—) - e f ^^5k(—). (5.33) 

!■ Us J J ° Us Jui=eo ^ Us 


We treat the interior term in (5.33), and place the boundary terms in the boundary 


contribution, /3 a ; which will be treated in the next subsection: 


9 

6 TV 

r'W'aj5a;<9i?( ) — 

Us 


-i^u^d^ Ve - Hr-hr vdR{ — ) 

r^ V lis "Ws Its 


(5.34) 


= d^^ l) + ( [5)^ 2) + ((5^3), 


(5.34 


(5.34 


^3/2 

1) = / / -+ e2 



^Uu,vd^ ( — ) ^ Ve\\u\\B, 

Us , 


2 ) = - / / —UujU^a^r^ ^ = -- j I —d^ (yl) ^ = - 


1 

2 Juj=9o Us 


Us 


2 < 


Us 

1 


u: 


'a. ( - 



„<5-l 


u 


2 Juj=9o Us 
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(5.34 


‘i) = ! I euojr^v^dn + e j j r^u^jvd^R ^ ^ 



For the final term in (5.34 


1 


3) , we have used: e sup i {R — Ro)\du}R{ — )| ^ C, which can 
_ _ [O,0o]4 


be estimated in a similar way as (5.30) - (5.31), to obtain: 


r^u^vdojR j ulr^ 



e\dujR{ — )\‘^v^r^] <Ve\\v\\l. 

Us ) 

(5.35) 


This establishes term (5.34) < \/e||u||^- 


now treat the second-order terms from (5.28): 


2 1 LJ 

.11 .) — 


ul-r^-^-e 


/u)=9o “s 


^VujLo ( — - 1 = (|5.36 


(5.36 


l) = e 2 


Us Ut 




1) (5.36 


2 ), 


/ We 

lu)=eo ^ Us 


(5.36) 




-^-11 r^-W 







2 ^Us 2 7^=0 Us' 


(5.36 


2 ) = e 2 


5—1 a I I 

r Vujdu, I 1 - e 

n I ^ 


Us J ^ JllJ=8o 

2 / 




lu}=eo 


5-1 Usu} 

r 

ui 


= € 


5-1 2 UsLO , 2 

^ Uuj—+e 

ui 


^v^vduiuj ' - e 





u 


lu)=eo 


5-1 Usu] 
r VujV^R. 

ui 


The middle term in (5.36 ,2) requires the following estimate on the profiles: 
sup j {R - Ro)\usuja}\‘^dR < C + esup j \ul^j‘^{r){R - Ro)dR 

<C + sup j \ul^ujir)\‘^rdr <C + 9Q^\\v\\jj3 <0o^e~^, (5.37) 


sup / I Us 


\R - Ro)dR < C + e'^ sup j \ulJ^{r){R - Ro)dR 

< C -b e0g ^||ug||vv2.4. 


r^'lulj^ + 


' ^eu.u.i 


(5.38) 
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t /■ 

Thus, we have term (5.36) < \/e| I'l’l |b+— / 

^ J U) 


The final second-order term from (5.28) is: 

Usuj 


5-1 2 1 f i5-l 2 1 2 f 

r V ^-—— / r v ^-—e / 

lu}=9o ^ Jui=0 J uj=9q 


5-1 ’Usu) 
r VujV^^. 

uj 


e I I r^~^^VRR 


— V- 


= (5.39 


1) (5.39 


(5.39 


~ V^ll^lls “ [ 

J U) 




= —e 

1 


VR 


2), (5.39) 

(1 -h 6)y/er^Vuj , 


Us 


+ 


Us 


+ r^+^dR{ — )vJ) 

Us ) 


ui=eQ 2 Us 


(5.39 


2) = e 


=. / / + (1 + 


-Fe / /. 


For the final term in (5.39 


(5.28): 


2), we use that e sup / r™'(i? — Ro)\usujr\‘^ < C. Thus, 

[O,0o] J _ 


(5.39) < ^/i||l;|| % — f -r^^^VR —. Finally, we have the low-order terms from |5.27| 
Joj=9o 2 Us _ 


e / ^rv r 2 r, / 1 2g3/2 / j.y r 2 O / 1 \ 

y/e -1- dR{rv) + r^vdR{ — ) -^/e-^- dR{rv) +r^vdR{ — ) 

^ \ Us Us Us J ^ \ Us Us Us 


+ C 


3/2 


i I Vcj Usuj . „ 

r vr{ - 1 - 2 

Us ui 


3/2 5-1 I 

e ' r Uui\ - 

Us ui 


Us 


ut 


2 J-1.. I »>- ) < C(e,«„)||„||| + iV||a||y 


(5.40) 


where C(e, 0o) —)• 0 as either 0o 0 or e —)• 0. Summarizing the results of this step: 
(ra + (ra <C'(e, 0 o)| 1^^111 +iV|klli + /^A, C{do,e)r.O{eo,Ve), 


(5.41) 


where /3a contains the boundary terms from (5.29), (5.34), (5.36), (5.39): 

.,2 


1 


Pa = -R 


— e 


2 Juj=9o 


luj=9o 


^ 1+^14 

Us ^ 2 


Usuj 5_i C' 

ui 2 



— e 

[ u^r^ ‘^dR{- 

J u=9q 

%j 

'UJ=Oo 

f „<5-l 2 1 

e 

f ...<5+1 ^ „,2 


~ 

/ r ^ —Vr. 

Jui=0 '^s 

2 

J u)=9o Us 


-) + - — 

Us 2 Ji^=Qq Us 

(5.42) 
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Step IV: Pressure Terms 


In this step, we apply our multiplier to the pressure terms from ( |1.25| — ri.27[ ), which 
immediately yield: 


Pcor^-^dR 


J-1 


2 

r V 


2 

.r V 


PRr^^^dU"-)= I I P.r^-^dnC—)- I I PRr^-^dU — ) 


2 

r V, 


u. 


Ju}=eo 

The interior term is estimated: 


ra„(—) + (S-1) I I ^er‘P 


U. ut 


Us 

(5.43) 


|< 5 - 1 || j j ^€r^p(^ 




Us Ut 


) I < 1^ —1| 


/p 2 


1/2 


tIIb- 


(5.44) 


We estimate the boundary term from (5.43) using the stress-free boundary condition in 


(1.32): 


Pp := [ t^-^POrC^) = -2e [ + 2e^P 

Jui=6o Ug Ju}=8n 




< -e 


luj=eo 


-ut + N 


U}=do “s 




luj=9o 


r^-^—+2e 

Us 


'u}=9o 


[ r^UujvdR{ — ) 

ui=9o 

(5.45) 


Step V: Boundary Terms 


We rewrite the Boundary contributions of the Navier Stokes terms, starting with (5.42): 

^2 r 


/3a = 

^ Juj=dQ Us 


1 2^-1 . f f 5-1 




2 Juj=0Q Us 


-uy-^ - e / u^r^-^dRi — ) + - I 

Juj=9o '^s 


2 Juj=9o Us 


— e 


/aj=6»o 


Usu} 5-1 £ / 

VVuj—r - — 

ui 2 


J-l 


a ;=0 

S-2^ 


1 

' Us 




,5+1 


1 


^ . o / 

" ./oj—, 


r —Vr 

u!=9q Us 


2 

,r V, 


2 Juj=eo 


,5-1 


-4 - e / ?■ UojdRi -) - e / 

J lj=6q j (jj=6o 


S-l UsLo ^ 2 

r vvu^-— V^ — 

^ J uj=0 ^ 


5-1 


^5+1 


2 Juj=0Q Us 


-Vr 


(5.46) 


2 Ji^=9q Us 


3/2 


5-lUuiV 

r - e 


lu}=eo 


Us 


Ilo=6o 


r^UujvdR{—) -e^ j 

Us Jui=0o 


5-1 Usuj 

r vvuj^r 

ui 


,5+1 


2 Juj=9q Us 


^ Juj=0 


7,2 

„5-l W 
Us ' 


(5.47) 
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In the equality yielding (5.46), the stress-free condition from (1.32) was used. Using the 
divergence free condition: — ^fkv = tvr ^ “ 2 -y/euU(^, we can rewrite 

two of the terms in (5.47): 


,5-1 


_ I - ul - - / ^ - vl = 

2 Us 2 Us 


I 


,5+1 


,5-1 


,5-1 


uj=9q “s 


2 Jui=9o 


-V + £2 


luj=8o 


-vUi^. (5.48) 


We also note that the in (5.47) is of a beneficial sign, and so we only keep treating 
the contributions. Summarizing, we have: 


/Sa = / r' 


5—1 U/^V 


— € 


/uj=0O 


Us 


[ r^Uajudni — ) - 

Jui=9o ‘^s 


3 
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,5-1 


Iu!=9o 

9)2 

5-1 W 


^_1 Usuj 6 

r vvo) 


2 2 


/ - - VUUJ - ^ / r- 

lu}=9o '^S 2 Jij—Q Us 


f U1+ 

luj=9o '^s 

(5.49) 


We now note that: / ev'^r^ ^ Using Young’s inequality we can 

J ijj=0n J J 


absorb all of the terms in (5.49) and (5.45) into either — j eu^r^ ^ or / ev^r' 


2„5-l 


/uj=9o 


terms except for the third term in (5.49), which we now estimate 


/llj=9o 


2 / 5—1 '^sui / ,5 Usui j 

-e r vvui^y = e r urv^^ =-e \ 

J ui=9q j ui=9q Ju]=9q 


udR{r^v^) 


ut 


= —e 


f u6r^ f t^uvr^^ + e f r^uvduiRi (5.50) 

J uj=9q ^s ^ uj=9(i '^s u ui=9o \ J 


= (5.50 1) + (5.50 2) + (5.50 3). 


First, we have: 


(5.50 


l)<e^ 


u 


'uj=9o 


'Ld=9o 


<eel 


1 

evlV <eel\\v\\l. 


For the second term, we use the divergence free condition vr = —y/e^v — ^Uuj, yielding: 


(5.50 2) = 6 


- ^3/2 


,5-1 


5-1_ '^suj 


/ui=9o 


< .3/2 


<< 


U 


'uj=9o 


-uv + 6 / r uu, 
J ui=9o 
1 
2 


u% 


<01=9(1 


W +6 


u 


< 01=9(1 


U,. 


<01=9(1 


ew + 


u,. 


u,. 


<01=9(1 
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+ [ eul + N{6)9oe f f ul < 660 f eul + e9o\\v\\%. 

J lo=6q J J J ijj=6q 


For the third term, we have: 


(5.50 


3) < '/~e\\duR ( — ) 1100 


ul,r^ I 9, 


.,2„<5 \ ~ <c a |L,||2 


We have used that 


I ( M loo 


1111 00 H“ 11 11 00111100 ^ C -\- 6||'Ug||j^4 ^6 ^ • 


The results of this step may be summarized as: 


nS-i 


\Pa\ + \Pp\ 


< - 


Jui=0q 

Step VI: Right-Hand Side 


(5.51) 


-ul + C{9o,e)\\v\\l, C{9o,e) ^ 0{9o,Ve). (5.52) 


[ [ fr^dR{—)= f f fr^ (^/e— + r^R{rv)—+r'^v^R{—)\ 

J J Us J J \ Us Us Us J 

<N{ 6 ) [ [ + S f f ev'^r^ + 6 [ [ r^\dR{rv)\‘^ , (5.53) 


e / / < N{ 6 )e 


< 7 V+^ + 5 


< N{ 6 ) / / eg‘^r‘^^^ + 5 


..2„<5 


ev r 


(5.54) 

(5.55) 


Placing the above steps together finishes the proof of Theorem 1.13 


6 Pressure Estimate 

Given P{lj,R) G L^(flAr), there is a corresponding scaled p{Lj,r) = P{uj,R), whose domain 
is n = (O,0o) X {Rq^Rq + \/fV). Abusing notation, denote the Euclidean counterpart 
to p{u),r) by p{x,y). 

Definition 6.1. Lq denotes the mean-zero subspace of L^: go G Tq iff f f qodxdy = 
f f qordrdio = 0. 
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Claim 6.2 (Mean-zero solvability of div). Denote the annular domain = (O,0o) 
{Ro,Ro + 9 q). For each go £ -^ 0 (^ 60 )’ t^iere exists a vector held vq G Hq{Q 0 ^) such that 
div ( Vo) = go ,l|vo||iTi(Oeg) < C'o|ko||L 2 (f^g^), where Co is independent of small 60 . 

Proof. The solvability of div : Hq — Lq is well known (see |Orlt98l pp. 26-28]). The 
important point for our analysis is that the constant Cq is independent of small Oq. This is 
guaranteed by |Galdill[ p. 162, estimate III.3.4], in which it is shown Co < ^ 

where R is the radius of a ball Bn C with respect to which is starlike. In our case, 
R ~ diam(D 0 (j). The claim is proven. □ 


Claim 6.3. For each g G L^(D 0 q), there exists a vector held v such that v = 0 on 
{R = Ro,Ro + 6 * 0 }, {w = 0 }, and < C||g||i 2 (ng^), where C is independent of 

small 60 . 

Proof. Similar to |Orlt98[ page 27] , dehne the mean-zero function qo = q—(^ f qrdrduj^ div(w), 

where w = ^60(]"^(r — 7?o)(^o — r + Ro)i^, 0^. By direct computation. 




dW{w)rdrduj = I, ||div(w)||i2 < \ I|w||hi < ^*0 


( 6 . 1 ) 


Denoting by vq the vector held guaranteed by Claim p?2\ for the function go, 
Ivoll^i ^ lko||L2 ^ lkllL2 + ( j j |g|rdrda;)||div(w)||i2 < ||g||L2||g||i26'o6'(7^ < 


The factor of 9 in the hnal inequality in (6.2) arises from Holder’s inequality: 


Igirdrda; <lkllL 2 


®0 


So rRo+So 


0 J Ro 


rdrdu) < ||g||L 2 ^'o- 


L2- 

( 6 . 2 ) 


(6.3) 


The desired vector held is now v = vq + ^ f f qrdrduj^w. Clearly v vanishes on the 
required components of the boundary, and we have: 


||v||j:^l < ||vo||h1 + 
The claim is proven. 



kk)l|w||/fi < ||g|U2 -h ||g||i26»o6»o^ = ||g||L2. ( 6 . 4 ) 


□ 


Claim 6 . 4 . There exists a vector held F(x, y) = {f{x, y),g{x, y)) such that div(F) = p{x, y), 
and: 


< C'IIpIIl 2 (o^^), (6.5) 
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where the constant C is independent of e,N, small 6q, and where F vanishes on the 
Dirichlet portions of the boundary {i? = Rq}, {oj = 0}, {ii = Rq + y/eN}. 


there exists a vector field, F^, such that div(Fk) = p on Fk(t<;, Rq + O^k) = Fk(a;, Rq + 
(/c + l)0o) = Fk(0, r) = 0, and ||Fk||/^i(Q)c) < C\\p\\i 2 (^^k^, where C does not depend on small 
00. Define F = J^Fk- Then ||F||^i(o) = < CY,M\LRn^) = C\\p\\l^, 

k k k 

and div(F) = p. Finally F satisfies the required boundary conditions. The claim is 

proven. □ 


Proof. Divide the domain into 0*^ = {0,0o} x {.Ro + ^^Oj .Ro +(fc + l)^o}. By Claim 


6.3 


The vector field F can be expressed in the polar coordinate basis ee and Cr, and as 
functions of uj,r, in which case div(F) = ^ + j + 'if r = Converting estimate (6.5) 

to polar coordinates reads: 

a2 j,2 


y/eN 


+ ^ +'^ + (j)^ +'ip'f) rdrdio < [ [ p‘^{u},r)rdrduj. (6.6) 
r r J J 


Dehne a(a;, R) = 4>{u}, r) and y/eb{ijj, R) = so R) = r), R) = 

4>rii^,r). Also, y/ebi^{Ld, R) = R) = 'ifriijJ-,^). Note that this scaling is the 


same as the Prandtl scaling. Scaling all of the terms in (6.6) yields: 


Oiv 


+ ^ + + -a|j + b]^ rdRdu < J J R(w, KfrdRdu, (6.7) 


and the scaled divergence equation: 


-h Ve- + bR = P. 

r r 


( 6 . 8 ) 


The admissible weights in (6.7) must be generalized to for 6 G [0,1]. The next claim 
shows this is possible as long as a small error is made in the scaled divergence equation 


( 6 . 8 ). 


Claim 6.5. Given P G L‘^{Qn) there exists a vector field Ai = (a, b) such that 

J J ^ r^dRduj ^ J J Kfr^dRdu, (6.9) 

where the constant is independent of N, e, and small 00- Ai vanishes on {R = Rq, Rq + N}, 
{cu = 0}, and satishes the scaled divergence equation: 


duj 

r 


+ V~e^+bR = P+( 
r V 


1-6 

2 


ebr 


-1 


( 6 . 10 ) 
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Proof. Given P, define P = Pra 2 . Applying the procedure culminating in estimate (6.7) 
to P gives a vector field (a, b) satisfying: 


-2 




and 


a“ + + ^ + + - 0 % + rdRdoo < 

\ T T C J J 


^ +Ve- + 6i? = P = Pr2-2. 


pV = 


PV^ (6.11) 


( 6 . 12 ) 


Define Ai = (a, b) = (ora 2 ^ ira 2 ). One readily computes the scaled divergence of Ai 


to check equation (6.10), as well as the desired estimates (6.9). It is also clear that A 


vanishes on {P = Rq, Pq + A}, {u = 0} because (a, b) vanishes on those components. The 
claim is proven. □ 


We now test against our equation against the multiplier {ar^,er^b) in several steps. 


Step I: Pressure Terms 

Applying the the multiplier {ar^,er^b) to the terms in equation (1.25 - 1.27) containing the 
pressure, P, yields: 


Pcoar^-^ + I I PRbr^ = - I I Pa^-^ + / Par^-^ - / / Pr%ii - 5 

U}=do 


J-l 


J- 


-P^r^ + 1(1-6) [ [ ^/^br^-^P + [ Par^-\ 

^ J J J u)=6o 


(6.13) 


We have used the relation (6.10). Using (6.9), the middle term above can be estimated: 

1 1 

, 1--5, 


Pr^-Web\<\l-5\ 


< 

rsj 


\i-6\eo 




p2^S 


2„(5-2 


2„<5-2 


eor' 


eb^f 


< 


(1-5)60 I I P^r\ (6.14) 


Step II: Laplacian Terms and Lower Order Terms 

In order to obtain the proper boundary cancellation, we use the representation of the 


Laplacian given in (4.3), (4.4). Applying our multiplier then yields: 


2e 


V ^ 3/2 ^ 

- Urr - Ur - fUojoj + -ffU - Vu) - VujR 


ar 


^V~ebP 
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+ 


2-v/e e 3-v/e e 

^VRR - VR --- 

1 ^ ry*^ ry^ 


ebr^ 


(6.15) 


We successively treat each term in (6.15), starting with the important, high order terms: 


URRT^a =11 r^URttR + I I 6r^ ^V^auR 


< 





1/2 



"^Ur] + J j 5r^ ^Ve-auR 


\ 1/2 


<^/i(^y j P{uj,Rfr^ 

< Ve\\P\\Ll5\\u\\A, 


1/2 



"^Ur] + j j dr^ ^V^auR 


\ 1/2 


(6.16) 


-2 


er^ ‘^auuioj = 2 / I er^ - 2e / ‘^auu 


ho=do 


< e 





1/2 





< (^\\P\\lI,5M\b - j 

J ijJ = 6 Q 


1/2 


- 2e 


“^auu! 


luj=9o 


(6.17) 


-2 



6r^ ^bvR 


(6.18) 


e^Vujujbr^ 2 


< e 


.5-2 _ / e\^br^-^ 



VujOujT —I c ( 7 (^L 

J uj= 9 q 
1/2 







1/2 




' U}=do 


<^M\B\\P\\Li,5-e^ [ ‘^Vu;b, 

J u)= 9 o 


(6.19) 


-e / / ^ = e / / r'’ - el '^URb 

J ui=9q 


„5-l„ 


„5-l„ 
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r^ul 


1/2 


< Ve 

< V~e\\u\\A\\P\\Ll,s - e 


eblr^-^ 


1/2 


- e / URhr 
J uj=9q 


S-l 


URbr' 


<5-1 


/ ui=9o 


( 6 . 20 ) 


The remaining terms can be estimated through Young’s inequality: 


—V^URar^ ^ + eaur^ ^ ‘^Voja + ev^jRar^ ^ + 2e‘^‘^bvRr^ ^ + Se^^'^bu^jr^ ^ "^bv 

< C(e, 0 o)||.f’|lL 2 (||ttm + ||i;||b) , (6-21) 

* ,0 

where C{9q, e) —)■ 0 as e, Oq —)■ 0. Summarizing the results from this step, 

4.3) X ar^ + (Equation 4.4) x br^ < C{6o, e)||P||j ^2 (||u|U + ||i^||b) 


(Equation 


— 2e / "^auaj — 

J u=9q 


*,(5 


‘^Vujb-e / URbr^ \ 


luj=Bo 


iui=eo 


( 6 . 22 ) 


where C{9o, e) —)• 0 as either argument —)• 0. 


Step III: Boundary Contributions 


We now treat the boundary contributions from the previous two steps. From (6.13) and 


(6.22), all of the boundary terms are: 


—e 


URbr^-^ - 


luj=9o 


[ “^v^jb — 2e / + f Par^ ^ = 0. (6.23) 

J u}=9q j u}=9q j uj=9o 


We have used the stress free boundary conditions from (1.32). 


Step IV: Profile Terms 


In this step, we give estimates on the terms from (1.25) - (1.26) which depend on the profiles 
Us, Vs- Precisely, we successively treat: 


1 1 \/c \/e 

-UsUuj + -UstoU + UsRV + VsUR H- VsU H- UgV 


ar 


+ 


11 2 1 

-UsVuj + -VsujU + VsVR + VsRV - ^UgU 

r r r Je 


ebr°. 


(6.24) 

(6.25) 
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Throughout the following estimates, we use the bounds on the prohles Us,Vs that were 

1/2 


proven in (4.19) - (4.26). 


UsU^ar^ ^<6*01^5 


u.,^uar^ ^ < lltt 



“'sa; oo 



. 2^6 





1/2 



1/2 

,2„5-2i < 




0o\\v\\B\\P\\Li,S, 


(6.26) 


1/2 


<^0 




\ 1/2 

j j UsRvar^ < 00 ^sup j - Rq)^ 


(6.27) 






1/2 


< 0o\\v\\b\\P\\l2,s, 


VsUrt a < llusri 


\/eVsUar^ ^ < \/e ' u^f 


(6.28) 




^<5„,2 
r Ufl 


.2„<5 


1/2 



1/2 



\ 1/2 

aV-2j <eo\\u\\A\\P\\Li,s, (6.29) 

1/2 



^Usvar^ ^ < 


UsVujbr^ ^ = lluslloo 


Vcujubr^ ^ < A/e 


^2^(5 



1/2 







, 1/2 


1/2 



<eoV~e\\v\\B\\P\\Li,5, (6.30) 

0o\\v\\b\\P\\l15, (6.31) 

1/2 

0o\\v\\B\\P\\Li,5, (6.32) 


,2„5-2 1 < 






^2+5 2 
^ ^R 


1/2 


1/2 






1/2 


VgVRbr^ < ^/e 

VsRvbr^ < ^oIIi^IIbII/^IIl^^, 

^/€Usubr^~^ < 

We summarize the results of this step: 

J J (Equation 




eb r 
\ 1/2 

<V~ee^o\\v\\B\\P\\Lls, (6.33) 

\ 1/2 

g^2^5-2j < \/e6»o||u||B||E’||i2_5, (6.34) 

(6.35) 

\ 1/2 

e6V-2 <02 ||i;||s||P||l2,5. (6.36) 


4.15) X ar^ + (Equation 


4.16) X ebr^ < (^(eo, e)||P||i2 (||u|U + ||ii||i?), 

*,o 


(6.37) 


where C{9q, e) —)• 0 as either 9q, e —>■ 0. 
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Remark 6.6. The choice of weight on the multiplier (ar*^, ebr^) is delicate in the sense that 
it is “critical” in several of the profile estimates given above. Specifically, in calculation 


(6.26), 11 Us 11 LOO is unable to absorb any factors of r, and so the weight after applying the 
multiplier (which in this case is UsU^ar^~^) must exactly match those of ||u||l and ||R||l2 ■ 


Step V: Right-Hand-Side 


Finally, we apply our multiplier {ar°,ebr°) to the right hand side of the system (1.25) - 


(|1.27|), which immediately yields: 

far^ + egr^b ^ j j ^ ^ ^ 

2„(5 , / / ,„2„(5 , r|iDll2 


„<5,r2 




< / //V+ / / 65 V+h||P||i.,,. 


(6.38) 


Putting the above steps together yields the desired Pressure estimate in Theorem 1.14 


7 Linearized Existence and Uniqueness for Navier-Stokes 
Remainders 


In this section, we prove Theorem 1.15 The full estimate for the linear problem, given in 


(1.50) is uniform in e, and small 9q. It was established in Corollary 2.8 that the spaces 


X and B endowed with their respective norms are Banach spaces, and so we can establish 
existence and uniqueness of the linear problem by applying Schaefer’s Fixed Point Theorem. 
We must apply the hxed point theorem for each fixed N, obtaining a solution u^,v^. We 


can subsequently send Ai —>■ oo as estimate (1.50) is uniform in N. 


For the following discussion, we fix an e and an < oo. Denote the Prandtl-layer version 
of the Stokes operator as 5*, so S^[u,v,P] = (/, 5 ) [u,v,P] satisfy the linear system: 


\/e e e 2 3/2 .It-, 7 

u - URR - Ur - kUuju + -^U - ' Vui + -Pui = /, 

V v/e e 2y/e e 1 „ 

- VRR - VR -H- 5 -Ua; + ^V + -Pr = g, 

^ 7* y*^ ^ 

Uuj + dR{rv) = 0 . 


(7.1) 

(7.2) 

(7.3) 


Claim 7.1. S-^ : L'^{nN) ^ U(Djv) is compact. 
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Proof. The usual Stokes solution operator, S~^, is compact on bounded domains from 
—)• L^. Suppose we take a sequence {fn,9n} £ ^{^n) which is uniformly bounded 

in the norm. Then define [/,g]„(a;,r) = [/n, ffnlCw, i?(r)), which are uniformly bounded 
in LP‘{PLi<[j) where N' = Rq + ^/eN. This implies [un,Vn\ = S{fn,gn) has a subsequence 
which converges in L^(f27v')- Define [un{uj, R), Vn{uj , R)] = [u{uj,r), ■^v{u;,r)]. Un,Vn solve 
the Stokes—* operator, so [un,Vn] = Sf^{fn,gn)- Moreover, by scaling, Un,Vn also have a 
convergent subsequence. Therefore 5"“^ is compact on 

□ 


We will need a version of Korn’s Inequality using the polar coordinate basis, for which we 
adapt the proof given in |Cia mi. 

Claim 7.2 (Lions’ Lemma). Let 17 be a bounded, open set with Lipschitz boundary. 
Suppose a distribution u G V'{U) has Vu = {^u,dru) G Then u G L‘^{U). 

Proof. Expressing ^ = costJUy — sincjUa, and Ur = cos ujUx + sincjUy, and the obvious 
inverse relationships, we have that {ux,Uy) G P[~^ [^u,dru') G P[~^. From here, the 

claim follows from Lions’ Lemma, found in [CialOt Thm 1.1]. □ 


Definition 7.3. Let e(u,v) be the symmetric gradient: 


e(u, v) = 


_ r 2 { r ^ 

hi^+Ur) Vr 


Denote the norm \\u,v\\e^ = \\u,v\\i 2 + \\e{u,v)\\E 2 + «:||%, ^||l 2 . 


(7.4) 


Claim 7.4 (Korn-type Inequality). For the solutions {u,v), we have the variant of Korn’s 
inequality: 


\\u,v\\h^u) ^\\u,v\\e^, (7.5) 

where the constant depends on the domain, U, but is independent of small k. 

Remark 7.5. As noted in |Cia m , these Korn-type inequalities do not make any restrictions 
on the behavior of {u, v) on the boundary dU. 

Proof. Using standard arguments, is a Banach space. We now show that Ei^{U) coincides 
with H^{U). Clearly, H^{U) C E^{U) continuously: ||'U,u||e« < C\\u,v\\e'^{u) where the 
constant is independent of small k. The reverse direction is delicate and requires a use of 
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Lions’ Lemma. We suppose ||u, < oo, so eij G =► Ve^ G Lf The components of 

V^('u,'u) can be expressed in terms of Ve^j via: 


- '^Lor '^Lur 
^ ? '^rr ? 5 

rp p 


j / '^(jj \ Uu) '^ojuj 


= 2 


du /- _ 

T 


'^ru) - 

r 


2 dr{ei 2 {u,v))+'^ 


'^LOr 

p2 


Thus, V^('U,u) G so by Lions’ Lemma, V(u,v} G L^, so (u,v) G The identity 
map i : is continuous and bijective, with bound independent of small k, and 

therefore by Banach’s inverse mapping the inverse map i~^ : ^ is also bounded. 

By observing that ||1 ||ek = I|1||hi = I|1||l 2([/), we see that the operator norm ||i“^||op 
independent of small k. 

□ 


Claim 7.6. 

2 


Z I UJ I Z I 

eu r + € -h Up^r + 


+ e^ — + ev'ir 

T ^ 


< 2 


(eu^ + e^u^)r+ j j 2e^ + 2ev\r + u\r + + / j 2euRVuj. (7.6) 


Proof. Expanding the definition of eij in (7.5), multiplying by 2 yields, and absorbing the 
^IIl 2 term to the left-hand-side gives: 




.U, 


2| |rt, u| \j_pi ^ “2^ J J ^ J J 2 Uj,r -|—— -|- Uj.r -|- 2 

Given our solutions {u,v), we define 

u{uj, r) = u{uj, R), v{u,r) = y/€v{uj, R). 

Scaling back to (w, R) coordinates yields the desired result: 


UyVijj . 


(7.7) 


(7.8) 


eu^r-b e-^-b n|jr-|- / / -b -b eu|jr 


,2 


< 2 


(en^ + eV)r+ / / 2e^ + 2eu|r + + 


2euRVuj. (7.9) 

□ 


Proof of Theorem\l.l!^ Consider the following map T{u,v) = (T^(m, u), T^(u, u)). 
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(7.10) 


T^{u, v) := f - \^UsUuj + -UsujU + UsRV + VsUR + —VsU + —UsV] + u, 


r 

2 1 


T^{u,v) := g - [-UsV^j + -VsujU + VsVr + Vsrv - -^Usu] + 


T is a bounded, affine map from Hi —)• Our solution is a fixed point of ^T, which 
is a compact map from Hi —)■ Hi. Let u^,v^ denote solutions to: 

n Ai7/^ 7/^ ^7,^ ^7,^ ^ 7/^ ^ f:3/2 A , 

o - A)u - URp. - —Uji - + 12 ^ “ 72 ^ ~ 


= A/ - A[-Usn^ + + UsRV^ + VgU^ + —VgU^ + —UsV^], 

r r r r 


. - . \ \ ^ A 2a/6 \ ^ A 

( 1 “-^)“- "^RR H 


A 


1 


j^2 ^ rp2 


e 

2 1 


= Xg - X[-Usvll + -VsujU^ + VsV^ + VsRV^ - :^UsU^\. 

r r r He 


(7.11) 


(7.12) 


We obtain bounds uniform in A in the following way. Select some 0 <Ao<<e<l. For 
Aq < A < 1 , the energy, positivity, and pressure estimates in the previous sections can be 
repeated to obtain uniform bounds (where the size will depend on Aq). For 0 < A < Aq << e. 


we must only perform an energy estimate by applying the multiplier {ru,erv) to (7.11 
7.12). For the upcoming calculation we drop the superscript on v^. 


-URR - ^UR - 2e^ ]ru= u^rV + 2e^ - 


UT. 


2e- 


UUn 


r Joj=eo r 
{-2vrr-—VR- e^)evr = 2e J J v^r + e — - e 


(1 - A) J J{u^ + v‘^)r + 

_ 3^3/2 I I ^ 

3 f /‘,3/2!f^_^3/2^^^<^ 

-ev^RU - eu^RV = 2e 


2e^v‘^ + eu^ 


7i2 \ 2 

.2W 


2\ I 


v^r 




URVoj - / tURV, 

J u)=9q 


PojU + PrVV = / Pu. 

J U) = 0q 


4^^ 


(7.13) 

(7.14) 

(7.15) 

(7.16) 
, (7.17) 

(7.18) 

(7.19) 
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The boundary contributions from (7.13, 7.14, 7.18, 7.19) cancel, using the same calculation 


as (4.11) and (4.29): 


/ Pu- euRV - / 2e— 

! u=9q J uj=9q Ju)=9q ^ 


— e 


vv,^ 


r 


= 0 . 


Ilj=9o 


(7.20) 


Summarizing the interior terms, and applying (7.6) 

,2 


UrV + 2e— + 2evRr + e^ — + 2euRVu} + (1 - + v^)r + 

r r 


> 2 


where 
J = 3e3/2 


eu^r + e—+URr+ f f + e^—+ eu^r 


+ 


(1-A) 


2e^u^ + 


r 

.2 I „.2 


+ J 


{u^ + u^)r + J, 
(7.21) 



Ua,V VuiU 


< Ve 


r 


2 \ 5 


r 


u: 


2 \ § 


On the right-hand-side, we have 


. „ , r 1 1 \/e -y/e , 

A/ - X[-UsUui + -UsujU + UsRV + VsUR ^ - VgU - UsV\ 


< X I I f^r + X I I u^r + ^ + u%r -h v^r, 


ur 


V 

r 


(7.22) 


(7.23) 


X 1 21. 

Xg - X[-UsVuj + -VsujU -h VgVR -h VsRV - l=UsU\ 

r r r Je 


^ A / / g^r + X I / ^ + u^r -h VrV -|- v^r, 


evr 


and therefore since 0 < A < Aq << e, we have 

{u^ + v‘^)r + j j -he|VeUp^r < j j (^f + eg'^'^r, 


(7.24) 


(7.25) 


uniformly in A. An application of Schaefer’s fixed point theorem then shows there exists 
a solution u^,v^ in the function space with the weak norm which we estimated (weak 


84 






















in the sense of weights and in e), and by linearity of our equation, applying the estimate 
above implies uniqueness of this solution. With existence and uniqueness in the weak space 
in hand, we can bootstrap: each is also an element of the space 

because Qj\f is a bounded domain. Because this solution is an element of X, B, we can apply 
the strong estimate |1.50 which is uniform in e, N, and small Bq. We can therefore send 
—)■ oo to obtain a global in R solution which also obeys estimate (1.50). This argument 


results in the proof of Theorem 1.15 


□ 


8 High Regularity Estimates 


In this section we prove Lemma 1.16 the high regularity estimate for the solution u, v to the 
problem (1.51 - 1.54). Notationally, we will keep the —M, where M > 0, as a large negative 
exponent for e, not taking care to rename different exponents as it is inconsequential to the 
estimate we are proving. 


Proof of Lemma We rescale via: 


do 


u{u, r) = u{6ouj, 6oR{r))-, v{uj, r) = y/ev{6ouj, 9oR{r)); P{u, r) = —P{6ouj, 6oR{r)). (8.1) 


The equation satisfied by the normalized profiles is: 

U Uu)ui . q2^ q - I ^0 1 1/ \ 

-Urr - t'o-^ - n()oVuj H-= — J := j(w,r), 

2 el _ - el ^ _ 

-Vrr - t'o-w + + -^v + Pr = —^9 := g{u},r). 


( 8 . 2 ) 

(8.3) 


Hf Estimates 

Using the standard regularity theory for Stokes equation, we can obtain estimates for u 
and V away from the corners of U. Let ^(w, r) denote a cutoff function which is supported 
near the corners of the domain in such a way that (w, r) G supp(x) ^ r — Rq < 1. Define 
Xi{oJ,r) = 1 — x(^) Then the equation for (?ii,i;i,Pi) := xi • {u,v,P) is given by: 


-Am - ^ + — = xif - 2Vxi • Vn - Axiu - ^XIljV + ^XIljP 

2. Iti 2 — 20 q _ 


~ ~ ^o)~2 ~ 




(8.4) 
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2 _ 1 _ - _ _ _ -2 
-Atii + -^Uiuj + -i^vi + Pir = X19 - 2Vxi • Vu - Axiv + dr{xi)P + -^Oodujixi)ui 


, /2 200 N 


(8.5) 


+ (l-0o^)^ + (^^)ui.. 

rpZ 

As (a;, r) £ supp Vxi U supp Axi, we have r — Rq < and so by the standard Stokes 
estimate (with inhomogeneous divergence, see Remark 8.1): 

||2Vxi • Vu - Axiu - \xiujV + \xiujP - (1 - Oq)^ - -— ^viuj\\l^ < e“^||/,5llL2, 


( 8 . 6 ) 

- 2Vy;i • Vv - Axiv + dr{xi)P + ^Ood^{xi)ui + (1 - 0o)^ + (“—< e“^|I/, 5 |Il 2 . 


The e ^ in estimate ( 8.6 


8.7|) arises from Axi = A (x(a;, 


(8.7) 

Thus using the 


standard Stokes estimate: 

II^i|Ia2 + I|1’i||r 2 + llA||j:^i < e 


P + \\9\\p) 


< 


L2 + \/i||5||L2) , (8.8) 


for some potentially large power M. We have used the calculation, according to the 
definition of / in ( 8 . 2 ) 

Il/lli 2 = / j prdrdu) = y/e j j prdRdu; = J j prdRdu = ele-^P\f\\\ 2 , 

(8.9) 

and analogously for g. Defining the corresponding profiles in Prandtl variables by inverting 
) above: ui( 0 ow, 6 'oA ='Ui(w, r), ui( 0 ow, 0 o-R) = ^ui(a;, r), Pi( 0 ow, 6 * 0 A = ^A(<^t), 


( 8.1 


we 


lave: 


ui,vi\\jj2<e ^||ui,ui||^2 < e ^\\LV~^9\\lI- 


-Ml 


( 8 . 10 ) 


Since M can be arbitrarily large, the quantity appearing on the left of (8.10) above is 
independent of e. We can also obtain weighted estimates by repeating the above analysis 
for the equation for wir 2"^2 ^ v^r‘ 2 ^ 2 : 

||,ij < e-"||rV2+^/2 11 ^,, (8 11) 

For the weighted estimate (8.11), we use that r < C on supp5^xi,/c > 1. This establishes 
the desired estimate for ui,vi. 
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Remark 8.1. We apply the Stokes estimate with inhomogeneous divergence because the 
weighted, cutoff vector field (■uira+i ^ 2 + 2 ) has divergence: 

Ouj , X r‘ 2^2 _ 5 1 _ 5_1 _ 15 / u 

— (r2^2ui) -ui + r2^2tij^^ + r2 2v^ = r2^2 I —|- 

r r \ r 


XlrV 


s 

7* 2 


1 

2 til. 


( 8 . 12 ) 


Therefore, a term e u||r 2 

'^*,5 

which is in turn controlled by e~^ 


appears on the right-hand side of the estimate (|8.8|) 
II/, y/€g\\j ^2 ^^ by our uniform energy estimates. 


^*, 2+5 


Estimates: 

Let y ;2 = xi^^, so 


(lo, r) G supp(x 2 ) ^ r — Ro < yfe and R — Rq < 1- (8.13) 


Define {u 2 ,V 2 ,P 2 ) = X 2 ■ {u,v,P)- By |OS93j . the Stokes problem has an estimate: 

||'W2,'il2||//3/2-I-||.P 2 ||j^i /2 < e ^ (^||/||l 2 + \/e||5lli,2^ ■ (8-14) 


We must relate the norm of U 2 ,V 2 to that of its scaled counterpart [rt 2 , U 2 ](a;, i?), 
given again by inverting transformation in equation (8.1). 


Claim 8.2. Define the transformation T : —>■ to be dehned in polar coordinates via 
'I'(a;, R) = (lo, r). Denote by u{uj, R) = u{uj, r) = if o T. Then ||m| {^ 3/2 = 1112 0 '^\\jj 3/2 < 
||lt||ji^3/2. Here the constant depends on the derivatives of T (which in turn depend on e). 


Proof. The transformation dl : (a;, R) —?• (a;, r) is bijective and has derivatives which are 
bounded above and below as a map from D C —?■ D C M^. Indeed, in the polar coordinate 
basis: 


VT(a;,R) 


^ U) 

R 

^2 

^ U) 

R 


'Ll 



From here it is easy to see that T is Bilipschitz, and | det d'l, | det T| ^ are bounded above 
and below, keeping in mind that the coordinate basis are functions of u. 


We may decompose ||ri||j:^3/2 = ||rt||Hi + II'^IIhi/ 2 . Using the definition of weak derivative 
and an approximation argument, it is easy to see the usual chain rule holds, namely 
Vu{x) = V'u^(T(x))DT(a:). Moreover, since the derivatives of ik are bounded above and 
below, we have by the change of variables formula: ||u||j:^i = ||uo THji^i < ||lZ||j|^i. 
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We must now treat the portion. As shown in |DPV12) . the norm of any 

function u is equivalent to the Gagliardo semi norm, which is defined as follows (where we 
calculate n + sp = 2 + ^2 = 3): 


\u]^ = 


\u{x) - u(y)|" 


Therefore, 

[u]^ = [u o = 

< II detDTiloo 



iQJn \x-y\ 
|h(^'(x)) - h(T(y))|2 


dxdy. 


Q Jfl 


|x — y\^ 

|tt(x') — M(y')P 


dxdy = 
dxdy' < 


\u{x') — tt(y')P 

nJn\^-^{x')-^-Ky'r 

u{x') - h(y')P 


detD'I'l ^dx'dy' 


„-il2 


\x' — 2/'P 


dx dy = [u] 


We have used that | detZ?T| is bounded above and below and that T is Bilipschitz. 


□ 


8.2 


we have: 11^2, U 21 1^3/2 ^ ^ ^ll^^2, ?^2 |Ih 3/2 ^ e ^ (ll/> VesIlL^) . We 


By using Claim 

can arbitrarily weight these norms due to (8.13). This concludes the proof of Lemma 1.16 


We are now able to control the high-regularity quantities appearing in 
equation (1.20): 


□ 

\z given in 


Proof of Theorem ] 1. Ill We first address the ur term in term (1.20). Using interpolated 
Holder and the uniform energy estimates from estimate (1.50): 


URllL^q < ||'Ui?||^2 \\ur\\\/ ^11/, ^£511^2 ||'Uk||^4 

* n-i-m ^ 1 I u xj Or^ * 


\i-e 


''*,<3 + 0: 


* , 2+5 


(8.15) 


We have used that a/q < 5/2 in order to apply the uniform energy estimates. Here 
9 = 9{5'), where 0 —> 1 as 5' —> 0. To estimate the term (8.15), we decompose 


u = ui -|- U 2 as in Lemma 1.16 




< II“1r|Il4 +||^t2R||L4 


*,2++ *'2+ + 
-Ml 


+ lF2ij||L4 


<IIwr||l 4 + 1^211^3/2 <e II/, \/e9llL2 . 

* 2 -I— ^ ^ ^, 2 "po 

’ 9 


(8.16) 

(8.17) 


In the second inequality of (8.16), we have used that R is order 1 on the support of U 2 
and so is equivalent to L^. In the first inequality in (8.17), we have used the 77^/^ 
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embedding in In the second inequality in (8.17), we used Lemma 1.16 to estimate 


U2\\h3/^- By observing ui\r=R^ = uiuj\r=Ro = 0, and ui\uj=o = uir\uj=o = 0, we use 

\\f,V~eg\\LK..- (8-18) 


estimate (2.13) for um to yield: 
\\uir\\l‘>‘ 




Inserting (8.17) into (8.15) and multiplying by 64: 


*, 2+6 
2 


*,2+5 


ei\\uR\\^2g <el ^4 ®)||/,+7^11^2 


(8.19) 


We can take S' small enough such that ? — M(1 — 0) > 0. The Uu and vr terms in 


(1.20) are treated in an identical manner, after observing the relevant boundary conditions 


are respected by the cutoff quantity: viR\r=RQ = = 0, ^ 111^=0 = uiij|^=o = 0, 

ni|r=i?o = uiaj\r=Ro = 0, and Mi|aj=o = uii^\i^=o = 0. We now treat v^j, which is slightly 


different from the previous terms. Recall from Theorem 1.17 that the parameter /3 > 0. 
Through interpolated Holder, we have: 


\\y/evuj\\^2q^^ < live' 


ew|li2 


WVeVujW]/ 


q 


< ll/,\/effllr2„ JIVeh 


^*,2+5 


\l-9 
'■^LO I I ^4 


( 8 . 20 ) 


Again we decompose v = vi + V 2 where V 2 is supported near the corner of the domain as 

~^\\LVeg\\L\ 


in Lemma 1.16 and so: 


2/3 — 2/3 + llv^^2tj||L4 2/3 ^ 


*,2+c5 

( 8 . 21 ) 


We cannot immediately apply Lemma 


2.16 


to the ||-v/euiij+4 term because vi = 
*,-M 

Xi(a;,R)u does not necessarily satisfy the stress-free boundary condition at {u = 0o}- 
However, ui^(0O) R) = Xii^o^ R)vuj{0o, R) + X1 lo{0o, R)v{9o, R), and so a trivial modification 
of the proof of Lemma |2.16| yields the required result. 

□ 


9 Nonlinear Existence and Uniqueness for Navier-Stokes Re¬ 
mainders 

We now apply contraction mapping on the space Z. For this section, call L the linear 


operator in the linearized problem appearing in equation (1.25). 
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Theorem 9.1. Suppose Lu,v = f{u,v),g{u,v). Select the 5' > 0 guaranteed by Theorem 


1.17 


, and let p = the Holder conjugate to g = 1 + 5'. Then for 6 G (0,1) sufficiently 


close to 1 , and for 27 + k < we have 
\\u,v\\z < C{Us,Vs) 


1_2_k II II 211 119 

1 + 62 2 '^\\u,v\\ze2\\u,v\\z 


(9.1) 


Thus, the solution operator to the nonlinear problem (1.25) - (1.32) maps the ball of radius 
2C{us,Vs) in Z to itself. 


Proof. We apply Theorem 1.17 after choosing the parameter ^ and subsequently <5 in 
the interval 1 — ^ < J < 1 , which yields < ||/, ^/eg\\l 2 . f,g are given by: 

^*,2+5 

f = f - ei (-UsUuj + -UsujU + UsRV + VgUR + —VsU + —Usv) , (9.2) 

\r r rr 


2/l_l_ _ _21_ 

g = g - ei -UsVoj + -VgcjU + VgVR + VgRV - -^UsU 

\r r r Je 


(9.3) 


where f,g are defined in (1.28) - (1.29). f,g are used to estimate the ||tt||x, ||i’| 


B 


components of ||u, v\\z according to the linear estimate (1.50), and the profile terms in (9.2) 


(|9.3|) are required to estimate the high-regularity components of ||u, u|| 2 , according to 

7 7 

The factor of e 4 accompanies these profile terms because ea was used in the 

We 


Theorem 


1.17 


definition of the norm || • jjz, while only a factor of e 4 was required in Theorem 
now proceed to estimate \ \f, ^/eg \\]^2 in terms of ||u, u||^. 

■^*,2 + 5 


1.17 


From Theorem 11.101 we have: 

-27- 


^ (^J j dRduj + e J j R.y+^dRduj'^ < ^-27-1^3/2-^ ^ gi-27-/c_ (94) 


Here we use that 27 -|- k < Next we have y/eR^’^, as defined in (1.30): 
e / / r\ulful < ||Up|||ooe [ f < C(up)e^"'"|||u,u||||, 


< e (^sup J(u^J^(R - Rq)^ J J 

J I r2+^(u^ij)^u2 < (^sup J{R-Ro){ul,ji) 




evR < ''|||u,u||||. 


|2 
1 00 


'"|||u,u||||. 


(9.5) 

(9.6) 

(9.7) 

(9.8) 
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u^ < e^ '"llltt, ?;||||, 
6^;2 < el 


n. V 




We now estimate \/ei?^’^, as defined in (1.31) 
' r\ulfvl<\\ul\\l, 


evl < e ^\\\u,v"'^ 


Z-, 


< e\\u\\l^\\vlJ\l 2 < e ^ '"\\\u,v\\\l 


"''lll^wllll, 


rH<fu^<H\\lo / / u^<e-^\\\u,v 


\z^ 


\ \ I I OO 

J j r'^^^vlp.fv'^ < e ^sup J (R - Ro){VpRf^ j j 


\u,v\\\z. 


(9.9) 

(9.10) 

(9.11) 

(9.12) 

(9.13) 

(9.14) 

(9.15) 


For (9.12) we use that ||ui ^^||^2 < e 2 and ||u|||cx> < e 2 ||tt|||, according to (2.10) 


For (9.15) we have used the bound: 



{vUf+ I 


esup J{R-Ro){vlRf < ^/esup J {v^r)'^ < ^ 

Summarizing the linear components of f,g, 

-^\Ru,VeRv\\L\ +Ve\\R"’^,VeR^’n\L\ < C(u„u,) + e5-i-||u,u|U. (9.16) 

*, 2 - 1-0 *, 2 - 1-0 


a-'y-9\ 


For the nonlinear terms we recall the dehnition of || • \\z in (1.20), where q = 1 + 5^ 
p = , and 0 < 2 < a < ^. We also recall the low regularity embeddings in Lemmas 

and 12.101 


2.9 


^ 27+1 


„<5,.2,.2 27+1 



U 


.2q 



y2pr5p 




u 


Z) 


(9.17) 


^27+1 


2+(5 2 2 ^ 27+1 
^ V U r, < e ^ 



y2p^p-1^5p 



^V+“ <e^||u,u|||. (9.18) 


Here the inequality holds because we have selected 5' > 0 so that ^ 


^ 27+2 


s 2 2 ^ 27+2 
r u V < e ' 



r2Sy^ 





(9.19) 
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Nonlinear terms in g: 


^27+1 


r^u\l < 


27 <6^-11 



Se'||«,r;|||, (9.20) 

<e'y\\u,v\\%, (9.21) 
(9.22) 


Combined with (9.16), we now have: 


\\f,Veg\\L\.^^C{us,Vs) + e2 I i^Wu^vWz+ e^^\ 


^lll- 


(9.23) 


We now provide estimates for the profile terms in (9.2) - (9.3) 


1 

£2 


1 

£2 


2 

£2 


e^+i 


gl+2 


1 

£2 


e'+i 


gi+i 




2 

£2 


e^+i 


r iigU^ ^ ll'^slloo^^ 




2+5 2-2 ^ Jii ||2 

S €2 ||t>sr||oo 


„-.2 5 




.-72 5 
Uj^r , 




< €2 NttJW 


< £2 { sup / uiRr\R - Rq) 


ev\\ 


2 .„2/ 



r^/2,72 < .ill,, ||2 

^ 11 11 00 


„2+5 2-2 < .1+ill,, ||2 

r Vg^R ^ f ^ N^sNoo 




// 


2+5-2 


r u 


R> 




„i5„-^2 


rW < NudlLe^ 




/ /rV. 


4 . 


Thus, (9.24) -|- ... -|- (9.34) < e 2 ''ll/,-v/e 5 '||f 2 • This concludes the proof. 

* * 1 I *, 2+5 


(9.24) 

(9.25) 

(9.26) 

(9.27) 

(9.28) 

(9.29) 

(9.30) 

(9.31) 

(9.32) 

(9.33) 

(9.34) 

□ 
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Corollary 9.2. For S sufficiently close to 1, the solution operator of the nonlinear equation 
is a contraction map on the space Z, satisfying: 


- U^,V^ - V^Wz <C{Us,Vs) 


€2 (||u^,u^||^ + ||u^ - u'^,V^ - 


v‘^\\z 


1 _ 7 _Lr 111 9 1 9 I I 

+ 62 2 — U ,V — v \\z 


(9.35) 


By applying the contraction mapping theorem, we have proven Theorem 1.18 and therefore 


the main result, Theorem 1.3 


Acknowledgements: The author thanks Yan Guo for many valuable discussions regard¬ 
ing this research. 
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